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Measure  the  characteristics  of  boundary  layer  turbu¬ 
lence  for  water  flowing  through  pipes  under  controlled 
laboratory  conditions.  Study  the  character  of  turbulence 
components  and  the  nature  of  water  flow  noise  by  harmonic 
and  numexical  analysis  of  exi^erimental  measurements. 
Compare  the  findings  with  existing  knowledge  of  airflow 
and  with  theoretical  analyses  and  computations  of  flow 
noise.  Derive  a  method  for  approximating  water  flow  noise 
values  for  a  range  of  conditions. 


RESULTS 


1.  The  water  pipe  flow  stucfy^gave  experiraentalXy  de¬ 
rived  values  of  pipe  wall  pressure  fluctuations,  radial  pipe 
deflections  at  one  point,  and  pressure -head  variation  along 
the  pipe,  for  turbulent  flow  rates  from  0.65  to  9.6  fps  with 
the  experimental  test  pipe  section  insulated  from  external 
noise  sources. 

2.  Test  pipe  material,  test  pipe  length,  and  location 
of  the  measurement  points  (hydro|^one  location,  etc. , )  did 
not  influence  the  measured  values. 

3.  The  broadband  fluctuating  pressure  at  the  pipe 
wall  depends  directly  on  the  dynamic  pressure  and  some¬ 
what  on  the  Reynolds  number  as  related  by  the  relation 
p/g  ^  0. 079/^** .  The  value  of  h  (the  velocity  ratio  of 
fluctuating  radial  stress  to  the  wall  shear  stress  in  the 
constant  stress  region)  averaged  0. 9  for  the  experimental 
data. 


4.  Based  upon  analysis  of  232  spectra  cujrves,  the 
shapes  of  the  curves  are  nearly  identical  and  are  of  the 
form  expected  from  theory.  Before  the  peak,  the  slope  of 


the  curves  is  about  1 .0,  and  after  the  peak,  in  the  equilib¬ 
rium  region,  the  slope  is  about  -1.66,  resembling  the 
Kolmogoroff- region  slope  of  the  kinetic  energy  spectrum. 

5.  The  position  of  the  peak  of  the  spectra  curves 
depends  upon  the  flow  conditions.  Based  upon  study  of  232 
spectra,  the  position  corresponds  to  a  Strouhal  number  of 
0,49  and  to  a  wave  number  of  150  feet*^. 

6,  Spectral  measurements  with  the  micrometer  set, 
with  very  thin  pipe  wall  thicknesses,  approximated  hydro¬ 
phone  measurement  spectra. 

7c  Numerical  analysis  of  photographed  waveforms 
produced  autocorrelations  and  power  spectra.  Spectra  from 
the  numerical  analysis  agreed  with  those  from  third-octave 
analysis.  The  Taylor  microscales  of  turbulence  are  of  the 
same  order  of  magnitude  obtained  for  air  flow  by  other 
experimenters, 

8,  The  dissipation  values  for  water  flow,  based 
upon  wall  shear  stress  measurements,  are  of  the  same 
order  of  magnitude  as  Laufer’s  values  for  air  flow. 


RECOMMENDATIONS 


1.  Apply  the  derived  method  for  approximating  the 
magnitude  of  broadband  fluctuations  to  experimental,  de¬ 
velopmental,  and  service  use. 

2.  For  research  purposes,  use  test  hydrophones 
which  resemble  those  used  in  actual  service  applications. 

3.  Continue  pipe-flow  "water  tunnel"  experiments  as 
an  efficient  method  for  studying  water  flow  noise  and  turbu¬ 
lence.  Study  higher  flow  rates  and  response  frequencies 


under  the  type  of  laboratory  conditions  described  herein. 
Study  application  of  experimental  findings  to  component  and 
system  development. 
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INTRODUCTION 


This  report  covers  a  continuation  of  basic  research 
initiated  by  C,  N.  Miller  of  the  Navy  Electronics  Labora¬ 
tory  to  study  water  flow  noise.  These  studies  have  in¬ 
cluded  experiments  with  streamlined  bodies  falling  through 
water «  to  show  that  measured  laminar  flow  noise  values 
would  be  less  than  turbulent  flow  values.  The  falling  body 
noise  experiments  were  conducted  with  19  laminar  test 
bodies,  at  Lake  Pend  Oreille,  Idaho,  and  were  completed 
in  December  1962.^* 

Experimental  studies  of  noise  produced  by  the  flow  of 
water  through  pipes  have  been  performed  at  Northwestern 
University  and  at  the  U.  S,  Navy  Experiment  Station, 
Annapolis,  Maryland,  *  *  However,  these  studies  were  not 
concerned  with  the  nature  of  flow  noise,  but  rather  were 
designed  to  show  the  sum  total  noise  in  piping  systems. 

The  pipe  flow  experiment  described  in  this  report 
was  devised  for  the  study  of  flow  noise  produced  by  forcing 
water  through  brass,  plastic,  and  glass  pipes  at  controlled 
rates,  under  laboratory  conditions.  Flow  rates  of  0,65  to 
9.6  fps  were  obtained  without  extraneous  noise,  and  accu¬ 
rate  flow  noise  measurements  were  made. 


’^Superscript  numbers  refer  to  the  list  of  references  at  the 
end  of  this  report, 

*  *This  activity  is  now  designated  the  U,  S,  Navy  Marine 
Engineering  Laboratory. 


WATER  PIPE  FLOW  NOISE 
EXPERIMENT 

The  objective  of  the  experiment  was  to  produce  water 
pipe  flow  noise  over  a  range  of  flow  rates  through  pipe  of 
different  materials,  under  controlled  conditions,  and  to 
collect,  for  further  processing  and  analysis,  three  basic 
measurements  related  to  the  nature  of  boundary  layer 
turbulence.  T^Mise  measurements  were:  1)  hydrophone 
pickups,  2)  pipe  wall  fluctuations,  and  3)  pressure 
differential  along  the  boundary  layer. 


INSTRUMENTATION 

The  experiments  were  conducted  in  a  concrete - 
floored  building  remote  from  industrial  and  traffic  noise 
and  vibration.  The  building  provided  space  for  a  65-foot- 
Ipng  straight  and  level  pipe  run,  and  had  separate  rooms  to 
help  isolate  the  test  pipe  section  from  the  sound  of  equip¬ 
ment  elsewhere  in  the  building. 

Fresh  water  was  stored  in  a  pressure  tank,  and 
entered  the  test  pipe  section  through  a  globe  valve  and  a 
30-foot  straight-line  assembly  of  soft  rubber  tubing  and 
special  fittings  (fig.  1).  The  tubing  was  used  to  eliminate 
conduction  of  noise  and  vibration  to  the  test  pipe  section. 

A  manometer,  calibrated  in  poimds  of  water  flow  per 
minute,  was  installed  across  a  3-inch  to  ^-incb  reducer 
between  sections  of  rubber  tubing,  to  function  as  a  venturi 
flow  gauge.  The  reducer  and  other  fittings  and  connections 
in  the  water  line  were  specially  designed  and  fitted  to 
minimize  flow  disturbances. 

Figure  2  shows  a  typical  test  setup.  The  water 
storage  tank  was  pressurized  through  rubber  tubing  from  a 
zero-  to  100-psig  controllable  external  air  supply  to  regu¬ 
late  the  water  flow  rate  through  the  test  pipe.  The  test 
pipe  section,  located  in  a  separate  room  with  measuring 
and  recording  equipment,  was  connected  to  input  and  output 
i^-inch  ID  soft  rubber  tubing  by  special  tapered  fittings  to 


Figure  1*  Input  section 
experimental  setup. 


of  water  pipe 


flow  noise 
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Figure  2,  Typical  water  flow  test  setup 


maintain  a  constant  inside  diameter.  The  output  tubing  led 
into  and  was  coiled  down  inside  a  collection  tank  installed 
on  Fairbanks  Morse  platform  scales.  A  viscometer  and  a 
float-type  flowmeter  were  inserted  in  the  rubber  outlet  line 
near  the  receiving  tank. 

The  test  pipe  section  was  of  several  configurations, 
comprising  3/8-inch  pipe  of  either  glass,  brass,  or  poly¬ 
vinyl  chloride  high-impact  plastic.  Test  sections  were 
supported  on  2  inches  of  rubberized  horsehair  acoustic 
padding  in  a  solidly  supported  wooden  trough.  The  pipe 
selected  for  the  experiments  had  an  inside  diameter  of 
0. 494  inch  and  a  wall  thickness  of  0. 09  inch.  The  test  pipe 
section  was  instrumented  with  either  of  two  types  of  test 
hydrophones  and  with  an  electronic  micrometer  to  measure 
radial  pipe  wall  fluctuations.  Plastic  and  brass  test  pipe 
sections  were  drilled  at  two  points  to  provide  differential- 
pressure  manometer  taps  to  measure  shear  pressure  drop. 

The  hydrophone  elements  were  installed  flush  with 
the  inner  walls  of  the  test  pipe  sections,  using  special 
adapters.  The  smaller-area  hydrophone  was  a  barium 
zirconium  titanate  bimorph  bender  0. 37  inch  long  and  0. 1 
inch  circumferentially.  The  other  was  a  radially  polarized 
barium  titanate  ceramic  cylinder  0. 45  inch  long  and  0. 494 
inch  ID,  moiuited  flush  inside  a  brass  adapter  fitting. 

Details  of  hydrophone  design  and  calibration  are  contained 
in  Appendix  A. 

The  electronic  micrometer  (Micrometer  Set,  Mutual 
Inductance,  AN/tTSM-43)  was  rigidly  mounted  to  10-inch- 
long  test  pipe  sections  which  had  been  thinned  in  that  area 
to  0. 003 -inch  wall  thickness.  The  micrometer  probe  was 
inductively  coupled  to  a  flat  metallic  surface  on  the  pipe  to 
measure  fluctuations  in  the  separating  distance  as  small  as 
ItT'^  inch.  (A  flat  brass  disc  was  cemented  to  the  plastic 
pipe  to  enable  inductive  coupling.  The  micrometer  was  not 
used  with  the  glass  pipe. ) 

The  hydrophone  was  connected  through  a  preamplifier, 
amplifier,  and  attenuator  to  a  Tektronix  512  oscilloscope. 


a  Magnecord  P-60  single -channel  tape  recorder,  and  a 
voltmeter,  A  Hewlett  Packard  20  c/s  to  5  kc/s  oscillator 
and  a  General  Radio  random  noise  generator  were  connected 
to  the  recorder  and  oscilloscope  through  an  attenuator  and 
the  preamplifier  channel.  The  mutual  inductance  microm¬ 
eter  set  was  connected  directly  to  the  oscilloscope,  re¬ 
corder,  and  voltmeter. 

Arrangement  of  the  readout  and  recording  equipment, 
shown  in  figure  3,  allowed  simultaneous  visual  display  and 
recording  of  alternate  hydrophone,  micrometer,  time- 
reference,  and  random  noise  signals. 


Figure  3,  Instrumentation  setup 


PROCEDURE 


The  storage  tank  was  filled  with  fresh  water  and 
pressurized  from  a  constant  compressed  air  supply  at  a 
predetermined  regulator  setting  corresponding  to  a  given 
water  flow  rate.  Zero-flow  hydrophone -circuit  ambient 
noise,  oscillator,  and  white -noise  generator  calibration 
reference  signals  were  alternately  tape-recorded  (the 
recorded  zero-flow  ambient  noise  was  no  greater  than  the 
theoretical  circuit  noise  value  computed  for  the  hydrophone - 
amplifier  circuit). 

The  water  collection  tank  drain  was  closed  and  the 
water  storage  tank  outlet  globe  valve  was  opened.  As  water 
flowed  through  the  test  section,  visual  readings  of  water 
temperature,  viscosity,  flow  rate,  and  shear  pressure 
drop  were  recorded  on  data  sheets.  Micrometer,  hydro¬ 
phone,  oscillator,  and  random-noise  generator  signals 
were  alternately  tape  recorded.  For  each  configuration  of 
the  test  pipe  section  (table  1)  test  runs  were  conducted  over 
a  range  of  flow  rates,  and  measurements  and  reference 
signals  were  recorded  for  each  run. 


TABLE  1.  WATER  FLOW  TEST  CONFIGURATIONS 


Pipe  Material 

Pipe  Length 

Instrumentation 

PVC  Hastic 

10  feet 

Hydrophone,  Micrometer, 

Manometer 

3  to  30  feet 

Manometer 

Brass 

10  feet 

Hydrophone,  Micrometer, 

Manometer 

3  to  30  feet 

Manometer 

Accurately  controlled  water  flow  from  2, 79  to  50 
pounds  per  minute  was  achieved^  giving  mean  average  flow 
rates  over  the  range  shown  in  table  2  (to  achieve  the  lower 
flow  rates,  the  outlet  tubing  was  constricted  between  the 
test  pipe  section  and  the  collection  tank). 

TABLE  2.  WATER  PLOW  THROUGH  0. 494-INCH  ID  HPE* 


Pounds  Water /Minute 

2.79 

10.0 

20.0 

30.0 

40.0 

50.0 

Flow,  Feet /Second 

0.56 

2.005 

4.01 

6.015 

8.02 

10. 025 

Flow,  Knots 

0.332 

1.187 

2.374 

3.561 

4.748 

5.935 

♦Based  on  fresh  water  at  60 *F,  with  a  density  of  62. 422  pounds  per  cubic 
foot. 


Water  flow  boundary  layer  turbulence  was  observed 
in  the  glass  pipe  over  the  I'ange  of  flow  rates,  using  a  lens 
system  and  a  bright  light  to  illuminate  suspended  particles. 


DATA  REDUCTION 

Raw  data  were  recorded  on  125,  000  feet  of  broadband 
audio  tape.  These  data  were  reduced  to  quantitative  form 
in  several  ways,  and  were  integrated  with  numerical  records 
from  the  corresponding  data  sheets  of  each  run,  for  analysis 
and  study  of  the  characteristics  of  the  flow  noise.  The  ex¬ 
periment  provided  necessary  data  to  derive  autocorrela¬ 
tions;  power  spectra,  broadband  magiiitudes,  and  spectrum 
levels  in  several  ways;  and  boundary  layer  wall  shear 
measures  of  energy  dissipation. 

A  Polaroid  camera  attached  to  the  oscilloscope 
photographed  the  hydrophone  and  micrometer  signal  patterns 
(fig.  4),  These  patterns  were  analyzed  by  Fourier  methods 
to  derive  power  spectra  and  autocorrelations. 

The  audio  tapes  were  played  back  through  a  Sanborn 
broadband  amplifier /recorder  to  give  graphic  analog  traces 
of  rectified  rms  hydrophone  and  micrometer  signal 
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amplitudes,  together  with  oscillator  and  random  noise  refer¬ 
ence  signal  traces.  The  broadband  recordings  were  also  sorted 
into  selected  third-octave  filter  channels  of  a  Sanborn 
Model  150  six-channel  recorder  and  an  NEL  .hird-octave 
digital  analyzer.  These  provided  analog  traces  of  rms 
values  against  an  accurate  time  scale,  and  digital  printout 
of  one -cycle  bandwidth  dB  values. 


WATER  FLOW  RATE  31  WATER  FLOW  RATE  36. 6 

POUNDS/MINUTE  'i  3  FPS)  POUNDS/MINUTE  (7. 35  FPS) 


5. 0  MILLISECOND/CENTIMETER  a  5  MIUISECOMD/CENTINeTER 


1. 67  MILLISECOND/CENTIMETER  a  5  MILLISECOND/CENTIMETER 


0. 3  MILLISECOND/CENTIMETER  a  5  MILLISECOND/CENTIMETER 


.Figure  4,  Photographed  osc illoacope  tracings  of 
broadband  audio  signals  from  cylindrical  hydro¬ 
phone  pickups  of  water  flow  noise  in  plastic  pipe* 
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DATA  ANALYSIS  AND 
PRESENTATION 

Results  of  the  water  pipe  flow  experiments  were  re¬ 
duced  to  several  graphic  forms  from  which  representative 
data  were  selected  for  correlation,  analysis,  and  presenta¬ 
tion.  Polaroid  photos,  broadband  rectified  rms  analog 
traces,  single -channel  rms  dB  analog  traces,  single-cycle 
bandwidth  rms  dB  values  in  digital  printouts,  and  numerical 
records  from  data  sheets  were  used  in  numerical  and 
harmonic  analyses.  Results  of  the  analyses  are  discussed 
briefly  and  presented  in  graphic  form.  Theory  and  methods 
used  in  the  analyses  are  presented  in  Appendix  B.  Con¬ 
clusions  and  recommendations  are  presented  in  a  following 
section. 

The  characteristics  of  turbulent  flow  are  generally 
defined  by  characteristic  lengths,  correlation  functions, 
characteristic  velocities,  and  by  spectrum  components. 

The  nature  of  flov/  noise,  in  this  case  the  measured  pressure 
fluctuations  at  the  pipe  wall,  can  also  be  related  in  terms 
of  the  same  types  of  variables.  The  one -point  hydrophone 
and  micrometer  measurements  furnished  broadband  and 
spectrum  levels,  and  the  broadband  magnitudes  are  shown 
to  be  related  to  the  velocities. 


AXJ  TOCORRELATIONS 

Autocorrelations  of  pressure  components  were 
measured  at  the  same  point  of  the  pipe  wall,  and  in  the 
same  radial  direction,  but  differ  in  their  time  separation. 
The  autocorrelation  is  a  periodic  function  of  time  lag,  and 
is  a  Fourier  transform  of  the  power  spectrum.  The 
Polaroid  photos  of  sound  pressure  waves  were  used  for 
numerical  evaluation  of  the  autocorrelations.  Results  are 
shown  in  figures  C  and  6 . 
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TIME,  MILLISECONDS 


Figure  5.  Autocorrelation  of  wall  pressure 
fluctuat ions,  calculated  from  oscilloscope  phot 
for  7.33  fps  flow  rate. 


ao  LO  2.0  3.0  40 

TiftlE,  MILLISECONDS 

Figure  6,  Autocorrelations  of  wall  pressure 
fluctuations,  calculated  from  oscilloscope  photos 
for  4,7  and  7,44  fps  flow  rates. 
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POWER  SPECTRA 


Spectxoim  levels  were  obtained  from  Sanborn  analog 
records,  TODA  digital  printouts,  and  by  numerical  analysis 
of  the  Polaroid  photos.  Figures  7  and  8  show  representa¬ 
tive  measured  spectra  for  the  two  types  of  hydrophones  and 
for  brass  and  plastic  pipe.  The  curves  show  presstme  level 
plotted  against  wave  number.  Approximately  232  spectrum 
curves  of  pressure  fluctuations  were  obtained  from  experi¬ 
mental  recordings.  The  shapes  of  these  curves  vary,  but 
have  fairly  constant  rising  and  falling  slopes  and  wave 
numbers  of  their  peaks.  The  average  rising  slope  is  1.0, 
and  the  average  falling  slope  is  1.G6:  this  is  nearly  equal 
to  the  Kolmogoroff  slope  of  -5/3. 

Figures  9  and  10  give  spectrum  levels  of  pressure 
calculated  from  the  photographs,  for  various  flow  velocities 
with  both  types  of  hydrophones.  These  curves  are  of  the 
same  shape  and  peak  at  about  the  same  wave  numbers  as 
the  curves  from  the  Sanborn  and  digital  analyzer  records. 
For  low  flow  rates,  the  curves  from  either  source  rise  and 
fall  less  steeply  than  on  curves  for  hi^er  flow  rates  (at 
the  low  flow  rates,  the  Reynolds  numbers  may  be  too  low 
for  the  inertial  subrange  —  in  which  the  flow  is  determined 
solely  by  the  dissipation  rate  —  to  occur). 

Figures  11  and  12  show  spectrum  curves  of  pipe  wall 
fluctuations  from  the  micrometer  measurements.  The 
rise  and  fall  rates  are  generally  less  than  for  hydrophone 
signal  curves.  This  may  result  from  pipe  wall  stiffness 
limiting  the  correlation  distances  to  small  values.  Sharper 
correlation  curves  tend  to  produce  larger  high-frequency 
components  in  the  system,  and  they  tend  to  flatten  the 
spectrum  curves. 

Spectrum  curves  for  the  experimental  data  peak  with 
a  mode  near  the  wave  number  of  150  feetT^.  Tlie  ordinates 
of  the  spectrum  curves  can  also  be  in  terms  of  the  nondi- 
mensional  Strouhal  number,  S  »  nr/u,  where  r  is  the  pipe 
radius,  u  is  the  mean  average  flow  velocity,  and  n  is  the 
frequency. 


OI«-CYCLE  BANDWIDTH  PRESSURE  LEVELS. DB 
REFERRED  TO  a  0002  DYNE/SQUARE  CENTIMETER 


WAVE  NUMBERS,  k  =27rn/u.  CONVERT  TO  STROUHAL  NUMBERS,  S  =  nr/u 
WHERE  n  IS  THE  FREQUENCY  IN  C/S 
r  IS  THE  PIPE  RADIUS  IN  FEET 
u  IS  THE  MEAN.  AVERAGE  FLOW  VELOCITY  IN  FPS 
BY  MULTIPLYING  THE  WAVE  NUMBER  BY  0. 3279  x  10‘2 


iO  100  1000  10,000 

WAVE  NUMBER,  FEET'^ 

Figure  7.  Cylindrical  hydrophone  sound  power 
spectrum  with  brass  pipe  flow. 
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WAVE  NUMBERS,  k  -  27rn/Q.  CONVERT  TO  STROUHAL  NUMBERS,  S  =  nr/Q 
WHERE  n  IS  THE  FREQUENCY  IN  C/S 
r  IS  THE  PIPE  RADIUS  IN  FEET 
u  IS  THE  MEAN.  AVERAGE  FLOW  VELOCITY  IN  FPS 
BY  MULTIPLYING  THE  WAVE  NUMBER  BY  0. 3279x10-2 


WAVE  NUMBER  FEET“^ 


Figure  8,  Bimorph  hydrophone  sound  power 
spectrum  with  plastic  pipe  flow* 
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ONE  -  CYCIE  BANDWIDTH  PRESSURE  tIVELS,  DB  REFERRED 
TO  a  0002  DYNE/SQUARE  CENTIMETER 


Figure  iJ.  Micrometer  signal  power  spectrum  for 
flow  rates  of  2,36,  4,52,  and  7,46  fps  through 
plastic  pipe. 


ONE-CYCLE  BANDWIDTH  PRESSURE  LEVEL,  DB  REFERRED  TO 
a  0002  DYNE/SQUARE  CENTIMETER 


Multiplying  the  wave  numbers  of  figures  7,  3,  9,  and 
10  by  0.  3279X1 CT®  will  give  the  equivalent  Strouhal  number 
for  the  0,494-inch  ID  experimental  pipe. 

The  maximum  of  the  spectrum  corresponds  to  the 
average  size  of  the  energy -containing  eddies,  and  the  eddy 
size  is  proportional  to  the  pipe  size.  The  spectral  maxima 
should  occur  at  the  same  Strouhal  number  for  pipes  of 
different  size,  since  both  eddy  sizes  and  Strouhal  numbers 
are  proportional  to  the  pipe  diameter.  When  the  spectra 
are  plotted  in  terms  of  the  Strouhal  number,  the  peak  should 
occur  at  5'  =  0. 49,  based  on  the  median  peak  value  of 
150  feet"^  in  wave  number.  The  ratio  of  the  pipe  diameter 
to  the  eddy  size  for  this  wave  number  is  6,18, 

The  spectrum  curves  are  plotted  in  the  same  di¬ 
mensions  as  the  kinetic  energy  spectrum  plots  of  Laufer, 
Schubauer,  Taylor,  Klebanoff,  and  others.  The  pressure 
spectrum  is  of  the  same  dimensions  as  the  energy  spec¬ 
trum,  Figure  13  compares  a  Laufer  spectrum  curve  with 
one  from  this  study,  with  the  curves  plotted  against  Strouhal 
number.  The  curve  from  this  study  has  a  wave  number  of 
120  feet*"^  at  the  peak,  and  the  Laufer  curve  has  a  wave 
number  of  0, 2  centimeter  at  the  peak,  Laufer  used  a 
9. 7 2 -inch  ID  pipe,  and  his  wave  numbers  convert  to  Strouhal 
numbers  by  multiplying  by  1,965,  The  Laufer  curve  is  for 
the  radial -direction  kinetic  energy  spectrum  as  measured 
at  a  distance  of  0.691  radius  from  the  pipe  wall.  The 
curve  Laufer  measured  nearest  the  wall  is  flatter  than 
other  Laufer  curves  and  most  of  the  curves  from  the  water 
pipe  flow  study  (probe  measurements  are  not  as  accurate 
neaj:  the  wsill).  The  Laufer  curve  and  the  curve  from  this 
study  are  similar,  and  both  peak  at  Strouhal  number  0.393. 

The  pressure  spectrum  curves  from  this  study  re¬ 
semble  the  energy  spectrum  curves  of  contemporary  fluid 
mechanics.  The  slopes  of  rise  and  fall  about  the  maximum 
are  nearly  those  expected  from  the  theory  of  isotropic 
turbulence. 


NOTE:  STROUHAL  NUMBER,  S  =  rnlU,  WHERE  r  IS  THE  PIPE  RADIUS  IN  FECT 

U IS  THE  MEAN.  AVERAGE  FLOW  IN  FPS 
n  IS  THE  FREQUENCY  IN  C/S 

- TYPICAL  WATER  PIPE  FLOW  SOUND  PRESSURE.  MEASURED 

USING  CYLINDRICAL  HYDROPHONE  IN  PLASTIC  PIPE, 

AT  6. 28  FPS  FLOW  RATE 

- LAlfER  RADIAL  KINETIC  ENERGY  SPECTRUM  AT  0. 691  RADIUS 

FROM  WALL  OF  9. 72-INCH  ID  WIND  TUNNEL^ 


ao  ai  Lo  lao  mo 

STROUHAL  NUMBER 


Figure  13*  Comparison  of  Laufer'a  kinetic  energy 
spectrum  xaith  pressure  spectrum  of  water  pipe  flow 
measurements  0 


BROADBAND  MAGNITUDES 


Broadband  amplitude  values  were  obtained  directly 
from  rectified  rms  analog  tracings,  and  from  numerical 
integration  of  spectra  levels  and  numerical  analysis  of 
Polaroid  photos.  The  Sanborn  traces  of  hydrophone  and 
micrometer  signal  levels  were  compared  to  parallel  traces 
of  the  reference  oscillator  and  white -noise  generator  signal 
levels. 


In  figures  14  and  15,  ratios  of  broadband  magnitudes, 
measured  from  both  types  of  hydrophones,  to  measured 
dynamic  pressure  head,  p/g,  are  plotted  against  Reynolds 
number.  Theoretical  piq  values,  based  on  the  equation 
p/q  =  k^F  (where  F  is  Blasius*  relationship  between  the 
friction  coefficient  and  the  Reynolds  number  equal  to 
0.079/^**),  with /c  equal  to  values  of  0.9,  1.0,  and  1.44, 
are  shown  as  straight  lines.  The  value  of  k  is  the  ratio  of 
the  fluctuating  radial  stress  to  the  wall  shear  stress  in  the 
constant  stress  layer  (it  is  assumed  that  the  fluctuating 
radial  treasure  at  the  pipe  wall  may  be  nearly  equal  in 
magnitude  to  the  fluctuating  pressure  in  the  conitant  stress 
layer).  The  measured  k  value,  with  Ihe  hydroplione  mounted 
flush  with  the  inside  of  the  pipe  wall  aiid  responding  only  to 
radial  pressure  components,  can  be  expected  to  be  0.9. 

With  a  hydrophone  responding  to  the  mean  cf  the  fluctuating 
components,  the  k  value  should  be  1,44.. 

Effective  broadband  magnitudes  measured  by  hydro¬ 
phones  are  shown  in  figures  16  and  17,  plotted  against 
measured  dynamic  pressure  head,  q  -  where  p  is 

the  fluid  density  and  u  is  the  mean  fluid  velocity. 


^OADBAND  MAGNiTUDE/DYNAMiC  HEAD,  p/q 


NOTE: 

1.  POINTS  PLOTTED  REPRESEffT  VALUES  OF  p/q  COMPUTED  FROM 
MEASURED  DATA. 

^  STRAIGHT-LI^$  CURVES  FOR  VALUES  OF  k  SHOW  THEORETICAL 
p/q  PLOTS  BASED  UPON  THE  EQUATION  p/q  =  k^F  WHERE  F  IS 
BLASIUS'  RELATIONSHIP  BETWEEN  THE  FRICTION  COEFFICIENT 
f  AND  THE  REYNOLDS  NUMBER;  F  =  a  079/R-  ^5  k  IS  THE  VELOCITY 
RATIO  OF  FLUCTUATING  RADIAL  STRESS  TO  THE  WALL  SHEAR  STRESS 
IN  THE  CONSTANT  STRESS  LAYER. 


p  IS  THE  MEASURED  RMS  BROADBAND  VALUE 
q  •  l/2pU  2 

REYNOLDS  NUMBER,  R  -uDp/p 
WHERE  PIS  THE  FLUID  DENSITY 

u  IS  THE  /V£AN.  AVERAGE  FLUID  VELOCITY 
D  IS  THE  PIPE  INSIDE  DIAMETER 
^  IS  THE  FLUID  DYNAMIC  VISCOSITY 


IGDO  10,000  100,000  1,000,000 

REYNOLDS  NUMBER 


Figure  14,  Bimorph  hydrophone  broadband  magnitude 
to  dynamic  head  ratios,  plotted  against  Reynolds 
numbe r. 
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BROADBAND  MAGNiTUDE/DYNAMIC  HEAD,  p/q 


i 


NOTE: 

1.  POINTS  PLOTTED  REPRESENT  VALUES  OF  p/q  COMPUTED  FROM 
MEASURED  DATA. 

^  STRAIGHT-LI^£  CURVES  FOR  VALUES  OF  k  SHOW  THEORETICAL 
p/q  PLOTS  BASED  UPON  THE  EQUATION  p/q  « k^F  WHERE  F  IS 
BUSIUS'  RELATIONSHIP  BETWEEN  THE  FRICTION  COEFFICIENT 
f  AND  THE  REYNOLDS  NUMBER;  F  =  0, 079/R-  25.  k  |S  THE  VELOCITY 
RATI  0  OF  FLUCTUATING  RADIAL  STRESS  TO  THE  WALL  SHEAR  STRESS 
IN  THE  CONSTANT  STRESS  LAYER. 


p  IS  THE  A/lEASURED  RMS  BROADBAND  VALUE 

q  -  l/2pu2 

REYNOLDS  NUMBER.  R 

WHERE  pis  THE  FLUID  DENSITY 

u  IS  THE  MEAN.  AVERAGE  FLUID  VELOCIT/ 
D  IS  THE  PIPE  INSIDE  DIAMETER 
p  IS  THE  FLUID  DYNAMIC  VISCOSITY 
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Figure  15,  Cylindrical  hydrophor^e  broadband 
magnitude  to  dynamic  head  ratios^  plotted  against 
Reynolds  number. 


FlUCTUATiNG  PRESSURE  ^OAOBAND  MAGNITUDES,  DYNES  PER  SQUARE  CENTIMETER 


FLUCTUATING  PRESSURE  BROADBAND  MAGNITUDES,  DYNES  PER  SQUARE  CENTIMETER 
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DYNAMIC  PRESSURE  HEAD,  q  =  l/2pD‘ 

WHERE  pis  THE  FLUID  DENSITY 

U  IS  THE  MEAN,  AVERAGE  FLUID  VELOCITY 
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DYNAMIC  PRESSURE  HEAD,  DYNES  PER  SQUARE  CENTIfleER 


Figure  17,  Measured  effective  hydrophone  broad¬ 
band  magnitudes  of  pipe  wall  fluctuating  pressure 
plotted  against  dynamic  pressure. 


WALL  SHEAR  STRESS/ DYNAMIC  HEAD  RAHO 


The  wall  shear  stress /dynamic  head  dimensionless 
ratio  can  be  used  to  evaluate  the  fluctuating  pressure  to 
(tynamic  head  ratio,  piq.  The  wall  stress /inertia  head 
ratio,  t/s,  is  nearly  constant  over  the  range  of  water  pipe 
flow  experiments,  but  is  a  function  of  the  Reynolds  number. 
Wall  shear  values  were  measured  with  a  differential - 
pressure  manometer  over  lengths  of  test  pipe  during  the 
water  flow  experiments,  and  were  recorded  on  data  sheets. 
Figures  18  and  19  show  the  wall  shear  stress /dynamic 
head  ratio  plotted  against  dynamic  head  and  against  Reynolds 
number.  The  variation  ith  Reynolds  number  resembles 
the  Blasius  friction  equation  discussed  imder  theory  in 
.  Appendix  B. 


PIPE  WALL  DISPLACEMENT 

The  mutual  inductance  micrometer  obtained  mea¬ 
surable  fluctuations  from  plastic  and  brass  pipe  after  the 
pipe  wall  thickness  was  reduced  to  0. 003  inch  near  the 
probe:  measurable  displacement  was  dependent  upon  wall 
thickness.  Radial  displacement  above  quiescence  of  plastic 
pipe  with  0. 003-inch  thick  wall  is  shown  plotted  against 
broadband  wall  pressure  in  figure  20.  Power  spectra  for 
radial  displacement  of  plastic  and  brass  walls  are  shown  in 
figures  11  and  12,  plotted  against  frequency. 


WALL  SHEAR  STRESS.  T=4Pr/L2 
Wh)ERE2lP  IS  THE  PRESSURE  DROP  OVER  THE  LENGTH  L 
r  IS  THE  PIPE  RADIUS 

L  IS  THE  LENGTH  OVER  WHICH  4P  WAS  MEASURED 
DYNAMIC  PRESSURE  HEAD,  q  :  l/ZpQ^ 

WHERE  pis  THE  FLUID  DENSITY 

u  IS  THE  MEAN,  AVERAGE  FLOW  VELOCITY 


DYNAMIC  PRESSURE  HEAD.  DYNES  PER  SQUARE  CENTIMETER 

Figure  18,  Wall  shear  strass  to  dynamic  pressure 
ratios,  plotted  against  dynamic  pressure  head. 


WAU  SHEAR  STRESS/DYNAMIC  PRESSURE  HEAD,  T/q 


WALL  SHEAR  STRESS,  r=4Pr/L2 
WHERE  2&P  IS  THE  PRESSURE  DROP  OVER  LENGTH  L 
r  iSTHE  PIPE  RADIUS 

L  IS  THE  LENGTH  OVER  WHICH  iP  WAS  A/IEASURED 
DYNAMIC  PRESSURE  HEAD,  q  =  l/Epu^,  AND 
REYNOLDS  NUMBER,  R  =  DD  piitl. 

WI€RE  p  IS  THE  FLUID  DENSITY 

u  IS  THE  MEAN,  AVERAGE  FLUID  VELOCITY 
D  IS  THE  PIPE  INSIDE  DIAMETER 
p  IS  TL€  FLUID  DYNAMIC  VISCOSITY 


Figure  19,  Wall  shear  stress  to  dynamic  pressure 
head  ratiost  plotted  against  Reynolds  number. 
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WALL  PRESSURE  FLUCTUATIONS,  DYNES  PER  SQUARE  CENTUVETER 

Figure  20,  Micrometer  measurements  of  0 , 003- inch- 
thick  plastic  pipe  wall  displacement ,  plotted 
against  broadba?id  wall  pressure  fluctuations 
measured  with  cylindrical  hydrophone. 


37 


CONCLUSIONS 
AND  RECOMMENDATIONS 


This  study  of  water  noise  in  pipes  has  demonstrated 
the  feasibility  and  economy  of  "water  tunnel"  experimenta¬ 
tion  in  turbulent  flow  noise  research.  The  study  also 
demonstrated  the  applicability  of  the  theory  of  classical 
mechanics  to  water  flow,  and  enabled  the  derivation  of  an 
engineering  method  for  computing  values  of  water  flow 
noise  for  given  conditions. 

Conduct  of  water  flow  noise  research  under  laboratory 
conditions  enables  control  of  the  variables  which  is  possible 
in  no  other  way.  This  basic  research  establishes  criteria 
applicable  to  further  laboratory  studies,  research  under 
more  severe  conditions,  and  to  developmental  programs. 

The  "water  tunnel"  type  of  experimentation  should  be 
continued  to  study  higher  flow  rates  and  response  frequencies. 
Application  of  the  results  of  this  study  to  design  and  develop¬ 
ment  of  components  and  systems  for  research  and  service 
use  should  be  encouraged. 
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APPENDIX  A: 

EXPERIMENTAL  HYDROPHONES 


Since  hydrophone  responsiveness  depends  upon  size, 
material,  shape,  and  polarity  of  the  element,  the  accuracy 
of  the  measured  data  depends  upon  accurate  hydrophone 
calibration. 

The  cantiiever-bimorph  type  shown  in  figure  A-1 
would  not  be  expected  to  respond  equally  over  its  face  if  it 
were  calibrated  by  a  plane  wave;  its  best  response  area  is 
near  its  free  edge,  and  its  effective  length  may  be  less  than 
half  its  true  length.  The  bimorph  t3rpe  experimental  hydro¬ 
phones  were  calibrated  in  a  free  field  and  also  within  a  test 
pipe  section,  using  a  small  probe  reference  hydrophone. 
Calibration  resr'  '■s  are  shown  in  figure  A-2;  the  response 
measured  within  the  pipe  is  several  decibels  higher  than 
that  measured  in  a  free  field* 

The  radially  polarized  barium  titanate  cylindrical 
type  hydrophone,  figure  A-3,  resembles  those  used  in 
service  applications.  This  type  was  cstlibrated  by  plane 
waves,  with  results  as  shown  in  figures  A-4  through  A-7. 
Its  response  to  wall  pressure  fluctuations  can  be  expected 
to  be  similar  to  the  sum  of  the  signals  from  the  elemental 
areas,  added  statistically. 


figure  A-1,  Barium  zirconium  titanate  bimorph 
hydrophone  with  element  mounted  in  machined  slot 
flush  with  pipe  inside  wall* 


MEASURED  AT  SWEETWATER  CALIBRATION  STATION 

WATER  TEMPERATURE  RS^’C  OPEN  CIRCI^T  VOLTAC^  MEASI^EO  ATTtf 
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— - -  NO  PIPE 

IN  PIPE 


Figure  A-2,  Simcrph  hydrophone  receiving  reaponae 
calibration  in  a  free  field  and  inside  a  pipe. 


Figure  A-3 ,  F.adially  polarised  barium  iitanate 
cylindrical  hydrophone. 


MEASURED  AT  SWEETWATER  CALIBRATION  STATION 
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Figure  A^4*  Cylindrical  hydrophone  receiving 
response  calibration. 
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Figure  A-5,  Cylindrical  hydrophone  calibration 
directivity  patternf  5  kc/s. 
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3  AUGUSTIN 


SV^OWATER  CALI  BRATION  STATION 

WATER  TEMPERATURE  26X  TEST  DISTANCE  2  ACIERS 

0EPTH3:9»^S  <^=90* 

SCAi£  10  DB  PER  RADIAL  DIVISION  ROTATE  $ 


Figure  il-5*  Cylindrical  hydrophone  calibration 
directivity  pattern,  10  hc/e. 
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Figure  A-7,  Cylindrical  hydrophone  complex 
impedance  calibration,  aeries  measurement. 


REVERSE  SIDE  BLANK 


DISCLAIMER  NOTICE 


THIS  DOCUMENT  IS  BEST  QUALITY 
PRACTICABLE.  THE  COPY  FURNISHED 
TO  DTIC  CONTAINED  A  SIGNIFICANT 
NUMBER  OF  PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


Hqi 


I 


PQ7 


available  copy 


APPENDIX  B:  THEORY 


MECHANICS  OF  FLOW 

The  flow  of  water  throu^  a  pipe  can  be  explained  by 
classical  mechanics.  The  behavior  of  a  continuous  medium 
(which  the  water  is  assumed  to  be)  under  mechanical  and 
thermal  environmental  influences  is  controlled  by  the  laws 
of  conservation  of  mass,  conservation  of  energy,  and  non¬ 
decreasing  entropy,  and  by  Newton’s  second  law  of  motion 
in  the  form  used  in  fluid  mechanics  imder  the  name  of  the 
Navier  Stokes  equations.  The  first  law  is  the  equatb^n  of 
continuity.  The  next  two  are  the  first  and  second  lav/s  of 
thermodynamics.  The  thermod3mamic  and  the  mechanical 
behavior  are  coupled  through  the  equation  for  the  conserva¬ 
tion  of  energy,  and  through  the  equation  of  state  which 
relates  stress  with  strain.  The  thermal  and  mechanical 
problems  can  be  separated  when  the  rise  in  temperature 
and  the  volume  changes  are  small.  The  cross -coupling  by 
compressibility  is  usually  small  for  a  nearly  incompressible 
fluid  like  water.  For  the  pipe  flow  experiments,  the  ther¬ 
mal  and  mechanical  problems  can  be  separated.  Only  the 
mechanical  problem  is  investigated  at  length  herein. 

Energy  is  required  to  push  the  water  through  the  pipe. 
P^iction  opposes  the  motion,  and  the  water  does  leot  slip 
freely  at  the  pipe  wall.  Part  of  the  friction  energy  appears 
as  the  pressure  energy  of  flow  noissr  with  most  of  the 
energy  eventually  converted  to  heat.  If  the  water  could 
slide  through  the  pipe  without  losses,  it  would  be  called  an 
ideal  fluid.  An  ideal  fluid  does  not  have  dissipation  and 
cannot  support  shearing  stresses  in  motion.  All  moving 
real  fluids  show  dissipation  which  appears  entirely  as  shear 
energy,  since  there  can  be  no  permanent  volume  changes  in 
real  fluids.  The  most  common  theoretical  fluid  considered 
is  the  Newtonian  fluid  in  which  all  the  dissipation  is  ex¬ 
pressible  in  terms  of  strain  rates  alone.  This  fluid  is 
defined  by  a  simple  relationship  (below)  between  thu  shear 
stresses  and  strain  rates  tlirou^  the  fluid-characteristic 
coefficient  called  the  viscosity,  which  determines  the 
dissipation:  » 


3 


(1) 


where  r 


is  the  shearing  stress,  with 

is  the  coefficient  of  viscosity 
are  the  velocity  components 

is  the  strain  rate  in  shear,  with  i^J 

are  the  space  distances. 


In  a  moving  fluid,  the  stress  depends  on  both  the 
static  strains  and  the  strain  rates.  In  an  isotropic  elastic 
solid,  the  equation  of  state  shows  that  the  stress  is  related 
to  the  strains  through  two  elastic  constants.  For  a  fluid 
that  is  not  moving,  mean  stresses  are  related  to  mean 
strains  through  one  elastic  constant: 

=  B(£k~aT)  *  B{3  ,c-ar) 

where  is  the  mecin  stress 

is  the  mean  strain,  equal  to  i&/3 
B  is  the  bulk  modulus  of  elasticity 
A  is  the  dilation,  or  volume  change  per  umt  volume 
a  is  the  coefficient  of  volume  expansion 
T  is  the  temperature. 

Static  fluids  cannot  support  shearing  stresses.  The 
state  of  stress  for  a  moving  fluid  is  the  sum  of  the  mean 
stress,  or  hydrostatic  pressure,  and  that  resulting  from 
the  fluid  motion.  Because  of  viscosity,  the  moving  fluid 
also  has  shear  stresses  that  are  proportional  to  strain 
rates.  The  stress  in  a  Amoving  fluid  is  given  by  the  formula: 


(2) 
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(3) 


/  bu.  bU  . 
(  +  ^ 
\^J 


.  p 
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tj 


where  a 


ij 


6 


tj 


is  the  stress 

is  the  Kronecker  delta  (equal  to  1  when  i=J 
and  zero  when  i4J) 


ir-  the  viscosity 

is  the  hydrostatic  pressure 

are  the  velocities 

are  space  dimensions 

is  free  Kronecker  index  for  velocity 

is  free  Kronecker  index  for  distance. 


The  energy  dissipation  per  unit  volume  in  a  Newtoman 
fluid  due  to  mechanical  causes  is  given  by: 

is  the  energy  diss4  ation  per  unit  volume 
is  the  stress 

is  the  strain  rate  in  shear  ii^J) 

is  the  viscosity. 

The  problem  of  the  viscous  fluid  is  stated  in  the 
following  equations: 

Continuity: 


where  ip 
or 


iJ 


bp  hipu^^) 
h 


(5) 
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where  p  is  the  fluid  density 

is  free  Kronecker  index  for  velocity 
X,  is  free  fCronecker  index  for  distance. 


The  First  Law  of  Thermodynamics: 


^dt  ^  at  ‘  ~dt 


de  dA 

’'dt 


where  h  is  the  enthalpy  per  unit  mass 
s  IS  the  entropy  per  unit  mass 

t  is  time 

T  is  the  tempei  ature  in  degrees  absolute 

e  is  the  internal  energy  per  unit  mass 

A  is  the  dilation 

p  is  the  fluid  density 

is  the  heat  conducted  per  unit  volume 

fi  is  the  viscosity 

V.  ,  is  the  strain  rate 

C  is  the  specific  heat 

X  is  the  distance. 

%/ 


The  Second  Law  of  Thermodynamics: 


where  s  is  the  entropy  per  unit  mass 
x^  is  the  distance 

g  is  the  heat  conducted  per  unit  volume 
T  is  the  temperature. 
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Newton's  Second  Law  of  Motion: 


ha 


hx 


J 


where  F  is  the  external  force  per  unit  mass 

or.  .  is  the  stress 

X  is  the  distance 
u 

u.  is  the  velocity 
t  is  time. 


The  First  Law  is  the  energy  relation  for  the  motion. 

It  can  be  used  to  develop  the  equations  for  the  thermal 
boundary  layer  using  Fourier's  law  of  heat  conduction 
(which  equates  heat  conducted  in  a  direction  to  the  product 
of  the  space  derivative  of  temperature  in  that  direction)  and 
the  constant  of  conductivity.  The  Second  Law  of  Thermo¬ 
dynamics  states  that  the  net  change  in  entropy  in  a  region 
(the  entropy  generated  in  the  region  and  conducted  out  of 
the  region)  must  be  zero  or  greater.  Substitution  of  the 
complete  stress -strain  relationship,  equation  3,  into  the 
Second  Law  of  Motion  leads  to  the  fundamental  equation  of 
fluid  mechanics,  the  Navier  Stokes  equation: 


Ibx^ 


where  p  is  the  fluid  density 

L  is  a  velocity  component 
t  is  time 

F.  is  the  external  force  per  unit  mass 

Ij 

^  is  the  coefficient  of  viscosity 


IS  free  Kronecker  index  for  velocity 
is  free  Kronecker  index  for  distance 

K 

is  the  distance. 

SIMILARITY  CONDITIONS  OF  TURBULENCE 

O.  Reynolds  made  a  detailed  study  of  the  character 
of  flow  in  188S.  lVi«.Ay  oi  the  concepts  of  turbulence  are 
due  to  Reynolds  who  noted  that  for  similar  flows  about 
geometrically  similar  bodies  with  different  fluids,  different 
velocities,  and  different  densities,  the  forces  on  the  fluid 
particles  must  bear  a  fixed  ratio.  When  various  measurable 
quantities  and  dimensional  constants  are  completely  related, 
an  equation  can  be  written  in  terms  of  dimensionless  ratios 
to  describe  the  relationship.  In  order  for  two  flows  to  be 
similar,  the  ratios  must  be  the  same  for  the  flows.  For 
similar  flow  processes,  the  dimensionless  ratios  are  the 
ratios  between  the  various  typ^s  of  forces.  The  principal 
ratios  are  listed  below: 


inertia 

viscous 

on 

M 

Reynolds  number 

pressure 

viscous 

F 

T 

pressure /shear  stress 

pressure 

elastic 

B  po^ 

pressure  /  acoustic 

pressure 

inertia 

3  '•mXmm 

pr 

pressure /dynamic  head 

inertia 

elastic 

pv^ ,  Jl! 

B 

Mach  number 

inertia 

gravity 

y3 

Froude  number 

viscous 

inertia 

T 

pV^ 

shear/ dynamic  head 

where  p  is  the  fluid  density 

V  is  the  representative  velocity 
I  is  the  representative  dimension 
fi  is  the  viscosity 
p  is  the  pressiire 
T  is  the  shear  stress 
B  is  the  bull:  modulus  of  elasticity 
0  is  the  velocity  of  sound. 

In  the  pipe  flow  experiments^,  the  elastic  (or  com¬ 
pressibility)  effect  and  the  gravity  effect  are  negligible. 

The  Mach  and  Froude  numbers,  therefore,  are  of  little 
concern.  Reynolds'  similarity  concept  mentioned  above 
states  that  one  dimensionless  ratio  is  a  function  of  the 
others.  The  number  of  dimensionless  ratios  is  the  sum  of 
the  number  of  related  variables  and  dimensional  constants 
less  the  number  of  primary  dimensions  like  length,  time, 
mass,  and  temperature.  For  the  usual  flov/  conditions,  it 
can  be  shown  that  the  coefficients  pIpV^  and  rlpV^  are 
functions  of  the  Reynolds  number  alone.  The  dimensionless 
drag  and  lift  coefficients  of  hydrodynamic  bodies  are  a 
function  of  Reynolds  number  alone  for  wide  ranges  of 
application. 


TURBULENCE  FLUCTUATIONS 

Reynolds  noted  that,  as  the  rate  of  flow  is  increased, 
the  flow  becomes  unstable  at  a  certain  Reynolds  number. 
Depending  on  the  quietness  of  the  initial  flow,  turbulence 
appears  at  a  higher  Reynolds  number  when  instabilities  are 
not  damped  out,  but  tends  to  grow  to  the  state  where  irreg¬ 
ular  flow  components  are  present  in  all  directions.  The 
Reynolds  number  where  turbulence  first  occurs  is  about 
2300,  based  on  the  equation: 

puJD  tjJ) 

M  V 


(10) 
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where  u  is  the  mean  flow  rate 
D  is  the  pipe  ID 
fd  is  the  dynamic  viscosity 
1/  is  the  kinematic  viscosity,  nip 
p  is  the  fluid  density. 

Reynolds  introduced  a  method  of  averaging  that  is 
still  in  use.  Turbulent  components  oscillate  in  time  about 
the  mean  value  that  is  steady  in  time.  The  mean  value  is 
defined  by  the  equation: 

t+r 

u  ^  /  u{t)dt  (11) 

where  T  is  a  time  interval  large  enough  for  the  mean 
value  to  be  constant  in  time  t 

u{t)  is  the  velocity  component  at  a  point  in  the  fluid 

u  is  the  mean  velocity  at  that  point. 

The  mean  values  of  the  turbulent  fluctuating  velocity 
and  pressure  components  are  zero,  since  turbulent  com¬ 
ponents  are  the  difference  of  a  velocity  from  the  mean: 

u  =  u  +  u '  (12) 

where  u  is  the  velocity  at  a  point 

u  is  the  mean  velocity  at  that  point 
u '  is  the  fluctuating  component  at  that  point. 

The  mean  square  values  of  fluctuating  components 
are  not  zero,  in  general.  The  quadratic  mean  values  are 
obtained  by  the  equation: 

_  t+T 

=  u'‘dt  (13) 

t 
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where  u  ®  is  the  mean  square  value  of  velocity  at  a  point. 

Tlie  turbulent  components  are  assumed  to  be  continu¬ 
ous  throughout  the  flow,  with  the  fluctuations  being  small 
compered  to  the  flow  dimensions.  Thus,  the  fluctuating 
components  are  assumed  to  be  continuous  functions  of  the 
distant  c  and  time,  and  the  basic  equations  of  mechanics 
apply  c’uring  turbulence.  The  equations  are  written  in 
terms  ef  the  values  during  turbulence  such  as  u  =  IT  +  u 
The  moan  over  time  of  the  equations  is  taken  by  the  Reyn¬ 
olds'  method  of  averaging.  Terms  of  the  first  power  of 
the  fluctuating  components  are  zero.  Additional  quadratic 
terms  t.uch  as  appear.  Many  studies  have  been  made 
of  isotropic  turbulence.  The  mathematics  is  simple,  and 
many  turbulent  processes,  such  as  the  flow  at  tne  center  of 
the  pipe,  tend  to  be  isotropic.  In  isotropic  turbulence, 
there  is  no  correlation  between  components  in  diffei'ent 
directions  at  a  point,  and  the  components  are  equal: 

ui  ^  uw  ^  m  =  0  u®  =  (14) 


where  u,  u,  and  w  are  the  fluctuating  velocity  components 
in  three  orthogonal  directions. 

For  a  fluid  to  have  turbulent  shearing  stresses,  there 
must  be  correlations  between  the  different  components,  and 
such  quadratic  terms  as  Xiu,  "M,  and  W  are  not  zero.  This 
occurs  in  non-isotropic  boundary  layer  flow  where  the 
influence  of  viscosity  is  of  importance.  For  pipe  flow,  the 
correlation  betv/een  the  radial  fluctuating  component  and 
the  axial  fluctuating  component  is: 


u  ti  =0  at  pipe  center  and  at  wall. 

f  JC 


0.45  at  maximum. 


(15) 


59 


STRESSES  IN  TURBULENCE 


The  Navier  Stokes  equations  of  motion  for  a  point 
show  stress  components  due  to  turbulence  in  addition  to 
those  due  to  laminar  flow.  These  stresses  are  called  the 
Reynolds  stresses,  or  the  virtual  stresses  of  turbulence. 
They  are  the  quadratic  time  mean  of  the  fluctuating  com¬ 
ponents  multiplied  by  the  fluid  density,  and  have  both  nor¬ 
mal  and  tangential  components  in  a  boundary  layer, 
and  with  i-ij„  They  add  the  tensor  to  the  stress 
strain  relation,  equation  3,  to  form  the  turlaulent  stress - 
strain  relation.  The  total  stresses  for  an  incompressible 
turbulent  fluid  are  given  in  cylindrical  coordinates  by  the 
equations: 
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where  V  .  V  .  and  U  are  the  mean  flow  velocities  in  the 
r  X 

radial*  tangential*  and  longitudinal 
directions 


and  u,  are  the  fluctuating  turbulent  velocities 
in  the  radial*  tangential,  and  longi¬ 
tudinal  directions 


P 


is  the  hydrostatic  pressure 


^  is  the  dynamic  viscosity 

p  is  the  fluid  density. 


In  order  to  solve  the  Navier  Stokes  equations  and  the 
stress -strain  relationships  for  turbulent  flows,  there  must 
be  obtained  relations  between  the  mean  and  the  fluctuating 
components.  J.  Laufer,®  P,  S.  Klebanoff,®  G.  B.  Schu- 
bauer,*  and  others  at  the  Nationsd  Bureau  of  Standards; 
British  fluid  mechanicists  such  as  A.  A.  Townsend;®  and 
other  scientists,  particularly  from  Holland  such  as 
J.  O.  Hinze,®  have  for  many  years  been  conducting  both 
theoretical  and  experimental  studies  of  boundary  layers, 
especially  in  air.  They  have  reported  relations  between 
the  fluctuating  and  the  mean  components  in  terms  of  non- 
dimensional  ratios.  These  ratios  are  mainly  dependent  on 
the  geometry.  The  energy  balance  in  the  boundary  layers 
has  also  been  reported  upon  very  completely. 


STRUCTURE  OF  TURBULENT  LAYERS 

Turbulence  of  fluids  flowing  over  solids  always  leads 
to  dynamic  similarity  conditions.  The  ?aw  of  the  wall 
states  that  turbulent  flow  near  walls  depends  on  the  wall 
stress  and  the  kinematic  viscosity,  v,  A  boundary  layer  is 
composed  of  a  thin  region  called  the  laminar  sublayer,  a 
constant  stress  layer,  and  an  outer  layer.  The  laminar 
sublayer  contains  turbulent  components,  but  in  lesser  mag¬ 
nitudes,  and  its  nature  is  mainly  determined  by  its  vis¬ 
cosity  and  the  wall  stress.  In  the  constant  stress  layer, 
the  motion  is  almost  entirely  determined  by  the  shear 
stress  and  the  viscosity.  Most  turbulence  is  produced  in 


or  near  the  rBgion  of  constant  stress.  There  is  energy 
equilibrium.  The  turbulent  energy  dissii  ation  into  eddies 
and  eventually  into  heat  is  about  equal  to  ihe  local  produc¬ 
tion  of  turbulence  in  this  region.  When  the  first  law  of 
thermodynamics,  or  the  energy  equation,  is  written  in 
terms  of  the  turbulent  quantities  (the  mean  plus  the  fluctu¬ 
ations),  the  energy  equation  can  be  separated  into  mean 
flow  and  turbulent  equations.  These  equations,  like  the 
first  law,  are  energy  balances  equating  the  sum  of  all  tne 
components  at  a  point  to  zero.  Figure  B-1  shows  the 
energy  balance  of  mean  motion  in  a  boumiary  layer  along  a 
smooth  wall  with  zero  pressure  gradient.  There  is  inter¬ 
action  between  the  inner  and  outer  parts  of  the  boundary 
layei".  In  the  out^r  region,  mean  flow  kinetic  energy  is 
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Figure  1 •  Energy  balance  of  mean  motion  in  a 
boundary  layer  along  a  smooth  wall  with  zero 
pressure  gradient  (adapted  from  ref* 


retarded  by  the  work  of  the  turbulent  shear  stresses.  The 
gain  to  the  mean  motion  by  convection  is  about  equal  to  that 
lost  by  the  turbulent  shear  stress  work.  The  production  of 
turbulence  is  small  in  the  outer  region,  but  it  increases 
strongly  toward  the  wall.  In  the  inner  region,  the  produc¬ 
tion  of  turbulence  exceeds  that  due  to  convection,  with  the 
difference  beiiig  made  up  by  the  work  of  turbulent  shear 
stresses.  There  is  an  influx  of  energy  toward  the  wadi 
where  it  is  converted  mainly  into  turbulence  energy, 
originating  from  the  mean  motion,  but  is  in  part  diffused 
back  by  turbulence  into  the  outer  region.  Figure  B-2  shows 
the  turbulent  energy  balance  in  a  pipe  at  the  wall  region. 

The  energies  of  dissipation  into  heat,  viscous  transfer  of 
kint  tic  energy,  production  of  turbulent  energy,  and  the 
diffusion  of  kinetic  and  pressure  energies  are  shown. 
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pipe  flow. 


A 

ESSURE 

DIFUS 

ION 

— 

S 

••TUR  BU 
_ W 

If  NT  EN 
ODUCTIC 

ERGY 

)N 

% 

:  DiFi 

JSION 

DISSI 

NATION 

'  VISCC 

KIN 

_ 

)US  TRA 
ETICEM 

NSFER 

ERGY 

KIP4ETIC 

ENERGY 

1 

Energy  balance  in  wall  region  of 


Close  to  the  wall,  all  of  the  terms  are  of  importance;  dis¬ 
sipation  and  production  are  nearly  equsil,  but  opposite  in 
algebraic  sign,  as  are  terms  representing  diffusion  of 
kinetic  energy  and  pressure  energy.  Further  from  the  wall, 
the  main  terms  in  the  energy  balance  are  the  production  and 
the  dissipation  with  some  turbulent  diffusion.  In  the  outer 
region  of  energy,  balance  is  between  the  turbulent  diffusion 
energy  and  the  dissipation.  The  peak  values  of  the  terms 
occur  at  the  wall  end  of  the  constant  stress  region  at  dis¬ 
tances  corresponding  to  wall  Reynolds  numbers  less  than 
1 0  based  on  the  v/all  distance  smd  the  friction  velocity, 

JCgU/v.  The  friction  velocity  is  defined  by  the  equation 


T  = 

w 


ML 

L2 


where  t 

is  the 

V) 

P 

is  the 

is  the 

LP 

is  the 

r 

is  the 

wall  stress 
fluid  density 
friction  velocity 

pressure  drop  in  the  length  L 
pipe  radius. 


Wall  Reynolds  numbers  of  about  8  have  been  called 
the  end  of  the  laminar  sublayer.  The  constant  stress  layer 
of  energy  equilibrium  extends  from  the  end  of  the  laminar 
sublayer  to  Reynolds  numbers  of  about  100.  Figure  B-3 
summarizes  figuratively  what  happens  to  the  energy  flow  in 
a  boundary  layer,  indicating  the  extraction  of  the  energy  of 
mean  flov/  and  its  transfer  in  the  direction  of  the  wall.  The 
greatest  production  and  dissipation  are  shown  to  occur 
near  the  wall.  The  turbulent  diffusion  of  energy  toward  the 
outer  region  is  also  shown. 
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Figure  B-3 ,  Energy  flow  in  a  boundary  layer* 


TURBULENCE  COMPONENTS 

The  exi>erimental  studies  of  Laufer  and  others  indicate 
that  the  velocity  ratios  are  mainly  a  function  of  the  geom" 
etry  and  are  nearly  constant  in  the  constant  stress  region. 
Figure  B-4  shows  the  ratios  of  the  fluctuating  components 
to  the  maximum  mean  axial  velocity  plotted  against  the  wall 
distance  radius  ratio,  x^lr.  Figure  B-5  shows  the  ratios 
of  the  fluctuating  components  to  the  friction  velocity  plotted 
against  the  wall  Reynolds  number,  x^u^jv*  These  velocity 
ratios  are  nearly  constant  in  the  constant  stress  region. 


ms  FLUCTUATING  VELOCITY/tVlAXIMUM  WIEAN  VELOCITY 


VS/ALL  DISTANCE/RADI  US  X2/r 

WHERE  U  I S  THE  MAX  IMUM  MEAN  VELOC  !TY 

X2  I S  THE  D I  STANCE  FROM  THE  P I  PE  WALL 
r  IS  THE  PIPE  RADIUS 

Ur  IS  THE  FLUCTUATING  RADIAL  TURBULENT  VELOCITY 
IS  THE  FLUCTUATING  TANGENTIAL  TURBULENT  VELOCITY 
Ux  IS  THE  FLUCTUATING  LONGITUDINAL  TURBULENT  VELOCITY 


Figure  B’-4,  Relative  turbulent  intensities  in 
pipe  flow  (adapted  from  ref*  S), 
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RMS  FLUCTUATING  VELOCITY/ FkICTION  VELOCITY 


WALL  REYNOLDS  NUMBER,  Rw  =  u,X2/i; 

WHERE  P  IS  THE  FLUID  KINEMATIC  VISCOSITY 
u,  ISM  FRICTION  VELOCITY 
X 2  I S  Tl€  D I  STANCE  FROM  rf€  P IK  WALL 
Ur  IS  THE  FLUCTUATING  RADIAL  TURBULENT  VELOCITY 
IS  THE  FLUCTUATING  TANGENTIAL  TURBULENT  VELOCITY 
Ux  IS  THE  FLUCTUATING  LONGITUDINAL  TURBULENT  VELOCITY 


Figure  B"~5»  Turbulent  intensities  near  the  wall 
in  pipe  flow  (adapted  from  ref*  8)* 


The  distribution  of  the  turbulence  shear  stress  is  of 
interest.  The  shear  stress /friction  velocity  ratio, 
is  plotted  for  pipe  flow  near  the  pipe  wall  (fig,  B-6),  The 
correlation  coefficient  for  shear  components  can  be 
approximated  from  the  velocity  ratios  of  figure  B-4  and 
from  the  shear  stress /friction  velocity  ratios  of  figure  B-5. 
The  velocities  are  the  rms  fluctuations  for  the  components: 


m 


\imsmis%imiamwmcm.  -iTvirf/u* 


3 


u  u  u  u  u  u 

r  X  r  X  r  X  - 

u.  uu  u.  u.  rx 

r  X  V  V 


u 


« 


or 


u  u 
r  X 


u, 


u  .8x2.5 
r  jc 


.45 


(18) 


VmiREYMOLDS  Rw  =  u,X2/v 

Vm  STRESS/FRICTION  VELOCITY,  -Tpf/u,2 
S^ERE  P  JS 1HI  FLUID  KINEf.lATIC  VISCOSITY 
u,  IS  THE  FRICTION  VELOCITY 

X2  \sm  oisTAf#::£  fro^t  the  pipe  wall 

ux  IS  THE  FLUCTUATING  LONGITUDINAL  TURBULENT  VELOCITY 
Ur  IS  m  FLUCTt^TING  RADIAL  TURBULENT  VELOCITY 


Figure  S“6,  Fiat  rihut ion  of  turbulence  shear 

atreaa  near  the  wall  in  pipe  flow  ^adapted  from  ref,  8), 
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The  correlation  coefficient  for  the  shear  is  equal  to  approx¬ 
imately  0. 45  over  most  of  the  pipe.  The  ratio  of  twice  the 
sum  of  the  turbulent  kinetic  energies,  or  and  the 

shear  stress,  -u^u^,  is  nearly  constant  over  most  of  the 
pipe.  Figure  B-7  is  a  plot  of  this  ratio,  -u^u^/  against 
the  wall  distance /radius  ratio,  xjr. 


WALL  DISTANCE/RADIUS  :X2/r 

TURBULENCE  SHEAR  STRESS/TURBUUENCE  Klf^TIC  ENERGY  -  UfC^/qj^ 

WHERE  X2  IS  Tl€  DISTANCE  FROM  THE  PI  PE  WALL 
r  IS  THE  PIPE  RADIUS 
IS  THE  SHEAR  STRESS 

qj2  IS  TWICE  THE  SUM  OF  THE  TURBULENT  KIFeiC  ENERGIES.  Uj  AND  Uj 


Figure  B-7»  Ratios  between  turbulence  shear  stress  and 
turbulence  kinetic  energy  in  pipe  flow  (adapted  from  ref*  B), 
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Figure  B-8  shows  the  turbulent  kinetic  energy  and  the 
turbulent  shear  stress,  -u^u^,  plotted  against  the  wall 
distance / radius  ratio. 
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TURBULENCE  KINETIC  EKSRGY/FRICTION  VELOCITY,  qfluj 

TURBULENCE  SHEAR  STRESS/FRICTION  VELOCITY.  - 

WHERE  X2  \Sm  DISTANCE  FROM  Tl«  PI WALL 

qj2  IS  TWICE  THE  SUM  OFT^S  TURBULENT  KINETIC  ENERGICES,  u,ANDu, 
u.  IS  Tt«  FRICTION  VELOCITY  ^ 
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Figure  5-S.  Diet  ri but  ion  of  turbulence  kinetic  energy  and 
turbulence  shear  stress  in  pipe  flow  ('adapted  from  ref,  2), 
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TURBULENCE  SHEAR  STRESS/RiaiON  VELOCITY 


PRESSURE/ DYNAMIC  HEAD  RATIO 


The  experimental  test  hydrophones  measured  the 
fluctuating  component  of  the  normal  stress,  at  the 
cylindrical  pipe  wall:  Laufer  was  able  to 

measure  the  distribution  of  the  Uj,q^  accounting  for  the 
viscous  transfer  of  kinetic  energy  toward  the  wall  in  the 
energy  balanc-'j.  The  diffusion  of  the  pressure  energy 
toward  the  wall  is  the  closing  balance  in  the  pipe  flow 
energy  balsince  curves  of  figure  B-2,  The  normal  fluctu¬ 
ating  pressure  component  at  the  wall  can  be  expected  to 
approximate  that  in  the  constant  pressure  layer. 


Willmarth  measured  the  wall  pressure  fluctuations  at 
the  wall  in  a  4-inch  ID  wind  tunnel  to  speeds  of  Mach  0.  S,’ 
His  measurements  are  given  in  terms  of  the  ratio  between 
the  fluctuating  rms  pressure  and  the  dynamic  head,  » 
with  U  being  the  free -stream  velocity  outside  the  boundary 
layer.  The  ratio,  pl^pU^  ^  has  a  nearly  constant  value  of 
0, 0035  over  the  Jdach  number  range.  (Theoretical  justifi¬ 
cation  of  the  constant  ratio  is  not  attempted  by  Willmarth: 
the  normal  fluctuating  pressure  must  be  primarily  related 
to  the  wall  stress. )  There  is  a  viscous  transfer  of  kinetic 
energy  toward  the  wall  from  the  constant  stress  layer,  and 
the  pressure  energy  is  also  transferred  toward  the  wall. 

The  fluctuating  normal  pressure  in  the  radial  direction 
should  be  nearly  the  same  at  the  wall  and  at  the  constant 
stress  layer.  Ratios  of  the  same  general  magnitude  as 
that  of  Willmarth  are  obtained  from  Laufer's  velocity  ratio 
measurements. 

If  the  ratio  of  obtained  from  the  ratios 

between  the  rms  radial  velocity  fluctuation  and  the  maximum 
mean  flow  velocity,  taken  from  figure  B-4,  a  rms  wall 
pressure  to  dynamic  pressure  ratio  results: 


2u  ® 

ptq  =•  2  0. 04^  =  0.0032 


(19) 
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For  a  mean  flow  velocity  distribution  across  the  pipe  of 
l/7th  power,  the  ratio  of  the  average  mean  flov/  rate  to  the 
maximum  mean  flow  rate  is  0.817,  The  wall  pressure/ 
dynamic  head  ratio  based  on  the  average  flow  velocity  is 
piq  =  0.0032/0.817®  =  0.00478,  using  figure  B-4. 


The  wall  stress /dynamic  head  dimensionless  ratio 
can  also  be  used  to  evaluate  the  fluctuating  pressure  to  dy¬ 
namic  head  ratio,  p/g.  The  wall  stress /inertia  head  ratio, 
r/g,  is  nearly  a  constant  over  all  the  experimental  mea¬ 
surements,  but  is  a  function  of  the  Reynolds  number.  Fig¬ 
ure  B“5  shows  the  radiejL  velocity/ shear  velocity  ratios  that 
can  be  used  to  obtain  the  p/g  ratios: 


pIq 


since  =»  where  is  obtained  from  figure  B-5 

and  is  assumed  to  be  ~0. 9.  It  is  assumed  above  that  the 
normal  fluctuating  pressure  at  the  wall  is  determined  by 


the  kinetic  energy  transported  from  the  ''onstant  stress 


layer  without  much  change,  p  =*  ,  For  an  approximate 

ratio  of  wall  stress /dynamic  head  of  0.0075  for  a  Reynolds 


number  of  15,  000,  the  p/g  ratio  would  be  0,0075  x  0,9® 


=  0.00608. 

Innumerable  experimental  measurements  of  wall 
friction  in  pipe  flow  have  been  made  over  many  years;  the 
wall  stress  per  unit  area  is  given  in  terms  of  the  dynamic 
head  and  a  friction  coefficient:  =»  /pu  ^  /  2  where  u  is  the 

average  mean  velocity,  and  f  is  the  friction  coefficient,  or 
the  ratio  r^^/q. 

Many  experimenters  like  Nikuradse®  have  obtained 
relationships  with  the  friction  coefficient  and  Reyno]  ds 
number.  H,  Blasius®  suggested  a  simple  rvlationship  that 
is  accurate  for  Reynolds  numbers  from  3000  to  100,  000: 

/  =  0. 0797?'^/*  where  /  is  the  friction  coefficient  or  T^^/g 
ratio,  and  P  is  the  Reynolds  number  based  on  the  average 
velocity  and  the  pipe  diameter,  uD/p,  Tlius,  the  p/g  ratio 
can  be  estimated  at  the  wall  from  the  friction  coefficient 
and  the  velocity  ratio,  u  ®  /u  * ,  by  the  equation; 
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(21) 


p!l  =  /  • 


For  the  water  pipe  flow  experiment,  broadband 
fluctuating  pressures  were  recorded  on  all  the  test  ruris. 
This  noise,  recorded  on  magnetic  tapes,  was  analyzed  with 
Sanborn  equipment  to  obtain  broadband  rms  magnitudes. 
Broadband  rms  magnitudes  were  also  obtained  by  numerical 
integration  of  the  spectrum  levels.  Since  random  pressures 
combine  as  the  sum  of  their  squares,  the  numerical 
integration  was  performed  by  dividing  the  frequencies  into 
narrow  bands,  in  which  the  spectrum  level  was  relatively 
constant,  and  adding  up  the  sum  of  the  squares  by: 


pa  = 


A/  +•  p 
^  1  as 


'Af  +P  ® 

8  ^8a 


where  p  is  the  rms  sp^ectrum  level  or  per  cycle  bandwidth 
level 

f  is  the  bandwidth  over  which  the  spectrum  is 
nearly  constant 

P  is  the  rms  broadband  pressure. 

The  ratios  of  the  broadbana  rms  m  'ignitudes  to  the  dynamic 
head  are  presented,  for  magnitudes  obtained  by  the  various 
methods,  in  figures  16  and  17. 


TURBULENCE  DISSIPATION 

Much  of  the  knov/led^  of  turbulence  is  due  to  the 
experimentalists.  The  equations  in  terms  of  the  mean  and 
fluctuating  components  usually  cannot  be  solved  without 
relationships  between  them.  Various  experimental  ratios  of 
fluctuating  to  the  mean  components  as  a  function  of  the  wall 
distance  have  been  discussed.  The  energy  halances  must 
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be  made  with  measured  terms.  The  most  complete  and 
reliable  measurements  with  turbulent  itities  are  the 
pipe  measurements  of  Laufer,  His  mea^^  its  were 

made  in  air  at  speeds  of  3  and  30  meters  per  second  maxi¬ 
mum  mean  velocity,  corresponding  to  Reynolds  numbers  of 
50,  000  and  500,  000  based  on  the  maximum  average  veloc¬ 
ities  and  the  pipe  diameter  (9.  72  inches  ID).  By  a  hot  wire 
anemometer  technique,  he  was  able  to  measure  triple  ve¬ 
locity  correlations  such  as  and  to  obtain 

terms  for  the  energy  balance  equation.  He  measured  the 


velocity  derivatives 


and. 


by  assuming  local  isotropy,  he  was  able  to  e/aluate  turbu¬ 
lence  dissipation.  Figure  B-9  shows  the  turbulence  dissi¬ 
pation  in  pipe  flow  based  on  Laufer' s  measurements, 
expressed  in  nondimensional  numbers  as  c!P/  2u^®  where  c' 
is  the  dissipation,  1)  is  the  pipe  diameter,  and  is  the 
friction  velocity.  This  quantity  has  a  value  of  approximately 
22  at  1/10  radius  from  the  pipe  wall;  for  greater  distances, 
it  varies  hyperbolicaliy  with  distance  from  the  wall.  The 
average  value  appears  to  be  about  1 1  if  the  dissipation 
coefficient  remains  the  same  for  points  nearer  the  wall 
than  1/10  radius  (the  dissipation  was  not  measured  in  the 
constant  stress  layer  which  is  nearer  the  wall).  Laufer' s 
values  are  a  careful  evaluation  of  the  dissipation  for 
isotropy  of  the  dissipation  tensor,  .Bu  Jox.^x. .  The 
values  give  local  dissipation  for  points  in  the  pipe  outside 
the  region  near  the  wall. 


Laufer's  values  of  the  dissipation  function  can  be 
compared  with  the  shear  measurements  of  the  water  pipe 
flow  experiment.  If  one  assumes  that  the  viscous  dissipa¬ 
tion  is  small  compared  to  turbulence  dissipation  (as  can  be 
expected  at  higher  Reynolds  numbers),  the  shear  power  at 
the  wall  can  be  assumed  to  be  equal  to  the  turbulence 
dissipation  in  the  corresponding  volume: 


('irr^L  =  T  ^2urL 

A  10  dt 


-L 

JO  dt  Y® 


(23) 
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DISSIPATION.  LAUFER'S  NONDIMEMSIONAL  VALUE 


REYNOLDS  NUMKR,  R  =  UD/J/=  50,000 
WALLDISTANCE/RAOUIS  =  2x2/0 
TURBULENCE  DISSIPATION  =  €'D/2u,3 
WHERE  U  IS  THE  /IrtAXI/^U/W  MEAN  VELOCITY 
0  IS  THE  PIPE  INSIDE  DIAMETER 
P  IS  THE  FLUID  KINETIC  VISCOSITY 
X2  I S  THE  D I  STANCE  FROM  THE  P I  PE  WALL 
€'  ISTHEDISSIPATION  INL2/r3 
u,  IS  THE  FRICTION  VELOCITY 


where  f  ' 

A 


T 

W 

r 


is  the  average  power  dissipation  per  unit  mass 
is  the  wall  shear  stress 
is  the  pipe  radius 


L  is  the  pipe  length  from  the  velocity  expression  dL/dt 


p  is  the  mass  density  of  the  fluid  {p=mj  V) 

la  ij.3 

is  the  dimension  of  turbulence  dissipation  which 
is  a  kinematic  quantity  of  energy  per  unit  time 
per  unit  mass. 


The  equation  above  states  that  the  average  dissipation 
is  directly  proportional  to  the  wall  shear  power  as  expected. 
In  figure  B-10,  is  plotted  against  wall  shear  power  for 
water  flov/  at  7.  5  fps  through  a  0. 494-inch  ID  pipe,  with  a 
measured  wall  shear  stress  of  0.  396  pounds  per  square  foot 
and  a  mass  density  equal  to  62.422  Ib/ft®  divided  by 
32. 174  ft/sec^o  The  turbulence  dissipation  for  these  con¬ 
ditions  is  148/63  ft® /sec®.  The  nondimensional  dissipa¬ 
tion  coefficient  (defined  by  Laufer  as  €'Dl2u^^)  which 
corresponds  to  this  friction  velocity  and  dissipation  value 
is  32.  7,  This  is  of  the  order  of  magnitude  of  Laufer’s  values. 
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Figure  Approx  iir.at  e  turbulence  dissipation 

in  a  0,494- inch  ID  pipe. 
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SOUND  POWER 


A  portion  of  the  fluctuating  turbulence  power  can  be 
called  sound  power,  Sound  power  per  unit  area,  or  the 
intensity  of  a  sound  wave,  is  defined  as  the  average  rate  of 
flow  of  sound  energy.  For  the  turbulent  components,  the 
intensity  is  the  correlation  between  the  fluctuating  pressure 
and  the  velocity  component  in  the  same  direction.  The 
sound  intensity  in  the  radial  direction  is  the  term  or 
pu^® .  Laufer  measured  the  triple  velocity  correlation, 

over  the  pipe  and  gives  a  value  of  0.45  over  most 
of  the  wall  distances.  An  index  of  sound  power  can  be  ob¬ 
tained  from  the  ratio  of  the  radial  sound  intensity  to  the 
average  friction  power: 


p  u 

fp  r 

T  V 

w 


r 

T  U 

V) 


T  VuJ 
w  * 


(24) 


where  t  is  the  wall  friction  stress,  pu.  *  =  t 
w  ^  ^  w 

f  is  the  wall  friction  factor 
U  is  the  average  mean  velocity 

is  the  radial  velocity  component  rms. 


This  sound  power  ratio  was  measured  to  be  0.  026  as  an 
average  for  the  water  pipe  flow  test  runs. 


The  intensity  of  the  radiated  sound  due  to  turbulence 
was  discussed  by  M.  J.  Lighthill  in  his  classic  papers  in 
the  Proceedings  of  the  Royal  Society  of  London.  By 

approximating  for  the  covariances  in  Lighthill's  expression 
for  the  radiated  intensity,  Proudeman  has  obtained  an 
equation  for  the  radiated  poM^er  per  unit  volume  for  isotropic 
turbulence.^  Proudeman  considers  the  equation  repre¬ 
sentative  at  higher  Rejmolds  numbers: 
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(25) 


6/2 

_  38p(u»^) 

where  P  is  the  radiated  power  pjer  unit  volume 
€'  is  the  dissipation  per  unit  mass 
0  is  the  velocity  of  sound  in  the  medium 
u’  is  the  fluctuating  velocity  component. 


RADIATED  POWER 

The  portion  of  the  wall  shear  power  radiated  as  sound 
can  be  estimated  roughly  using  the  average  dissipation  and 
Proudeman's  equation.  The  fraction  of  the  power  radiated 
is: 
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.  rdL 
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(26) 


It  can  be  seen  that,  according  to  this  formula,  the  fraction 
of  power  radiated  into  sound  is  very  small,  except  for  high 
fluctuating  velocities.  For  a  fluctuating  velocity  of  10 
centimeters  per  second  (corresponding  to  a  water  flow  rate 
of  30  pounds  per  minute  through  the  test  pipe)  and  with  the 
velocity  of  sound  in  water  of  1.48x10®  centimeters  per 
second,  the  fraction  radiated  into  sound  is  28(10/1.48x10®)® 
=  0.53x10-20. 
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HARMONIC  ANALYSIS 


All  of  the  pressure  fluctuations  of  this  report  were 
measured  at  one  point  of  the  hydrophone  output  at  the  pipe 
wall.  The  photographs  of  both  the  hydrophone  and  the 
micrometer  outputs  show  the  random  nature  of  the  pressure 
fluctuations.  Tlxe  waveform  is  constantly  changing  its 
shape  with  time  and  never  exactly  repeats  itself  within 
finite  time.  The  infinite  number  of  waveforms  is  called  an 
ensemble.  Functions  of  this  type  are  described  by  statis¬ 
ticians  in  terms  of  probability  distributions  for  the  in¬ 
stantaneous  values  to  arrive  at  statistical  averages.  The 
averages  of  this  analysis  are  values  measured  by  tech¬ 
niques  based  on  harmonic  analysis.  Generalized  harmonic 
methods  are  suitable  for  analyzing  most  functions  such  as 
periodic,  transient  or  aperiodic,  and  random  functions.  It 
has  been  mentioned  earlier  that  turbulence  is  primarily 
described  by  a  representative  velocity,  a  representative 
length,  space  and  time  correlations,  and  power  spectra. 

The  present  analysis  principally  gives  the  broadband  rms 
magnitude  which  is  related  to  the  velocities,  and  the  power 
spectra.  Autocorrelations  are  also  analyzed.  They  play  a 
central  role  in  the  general  harmonic  theory.  The  power 
spectra  are  obtainable  without  the  autocorrelations.  Har¬ 
monic  analysis  is  based  on  the  assumption  by  the  Wiener 
theorem  that  autocorrelations  exist  for  the  functions  and 
that  a  common  autocorrelation  exists  for  all  the  members 
of  the  ensemble.  This  means  that  all  members  of  an 
ensemble  have  the  same  power  spectrum  regardless  of  the 
waveform.  Due  to  the  fluctuating  nature  of  turbulence  with 
zero  mean  values,  and  because  the  relationships  are  usually 
between  two  quantities  expressed  as  a  prodxict,  the  variables 
of  the  analysis  are  mostly  in  quadratic  terms  such  as  the 
mean  squared  values,  the  power  spectra,  and  the  correla¬ 
tions  (the  broadband  values  were  discussed  in  conjunction 
with  press’are/ dynamic  head  ratios). 


CORRELATIONS  DEFINED 


The  correlations  between  quantities  in  flow  have  been 
shown  to  be  tensors  by  de  Karman  and  Howarth,^®  Correla¬ 
tions  in  a  steady  flow  are  between  fluctuating  vector  com¬ 
ponents  integrated,  as  a  product  over  a  sufficiently  long 
time  period,  to  the  correlation  independent  of  time: 


‘  llmT- 


1  r 


where  If 


ij 


(x,  r)  is  the  correlation  between  vector  comp)onents 
i  and  J  displaced  by  a  vector  distance  r 


ix)  is  the  instantaneous  value  at  point  x  of  the  i 
component  of  the  fluctuating  vector  velocity 


r) 


is  the  instantaneous  value  of  the  fluctuating 
vector  J  component  at  the  point  displaced  r 
from  the  point  x 


T  is  the  interval  of  integration. 


Correlations  in  flow  have  been  defined  in  various 
manners.  Figure  B-11  shows  some  of  the  definition::.  In 
1937,  de  Karman  and  Howarth  defined  longitudinal  correla¬ 
tion  as  the  correlation  between  components  at  points  A  and 
B  in  the  same  direction  as  the  line  between  the  points. 
Lateral  correlation  is  defined  as  that  between  two  compo¬ 
nents  at  right  angles  to  the  direction  between  points  A  and  B. 
Taylor's  experimental  measurements  were  made  with  hot 
v/ires  which  responded  mainly  to  components  in  the  direction 
of  flow  in  the  wind  tunnel.^**  ^  Taylor's  longitudinal  cor¬ 
relation  is  between  velocity  components  in  the  direction  of 
the  main  flow  at  two  points  located  in  a  line  in  the  main 
flov/  direction.  Taylor's  lateral  correlation  is  between 
components  in  the  main  flow  direction  but  displaced  at  right 
angles  to  the  tunnel.  Harrison  at  David  Taylor  Model  Basin 
and  Von  Winkle  at  the  U.  S.  Navy  Underwater  Sound 
Laboratory  relate  pressure  values  at  the  wall  in  their  wind 
tunnels.  Pressure  correlations  are  called  longitudinal 


TYPE  CORRELATION 

EXP 

experimenter 

ERIMENTALISTS'  OEF 

CORRELATION 

INITIONS 

PRESSURE 

COA1PONENT 

EETVJEEN 

POINTS 

LONG  ITUD  INAL: 

COMPONENTS  ARE 
PARALLEL  TO  DIRECTION 
OF  DISPUCEMENTS 

ds  kArMAN  & 
hovvarth 

LONGITUDINAL 

LONGITUDINAL 

ATOB 

TAYLOR 

LONGITUDINAL 

LONGITUDINAL 

ATOB 

LATERAL: 

SINCE  PRESSURE 
COMPONENTS  ARE 
NORfMLTOTHE 

DIRECTION  OF 

THE  DisPL' ::em£nts 

da  KARMAn  & 
KOWARTH 

LATERAL 

..ORRIALTOWALL 

DTOE 

TAYLOR 

UTERAl 

LONGITUDINAL 

ATOC 

HARRISON 

LONGITUDINAL 

NORfMLTOWAa 

DTOE 

TRANSVERSE 

NmiALTOVm 

DTOF 

VON  WINKLE 

LONGITUDINAL 

KX)RMALTOWAa 

DTOF 

LATERAL 

NORfMLTOWAlL 

DTOF 

UlTIf^N 

AUTOCORRELATION 

NORMAL  TO  WALL 

GT06 

Figure  B-11,  Sketch  of  space  correlations  for  pipe  floio. 


81 


for  points  in  the  direction  of  tunnel  flow.  Harrison  calls 
correlations  between  wall  poirts  at  right  angles  to  the  tunnel 
transverse  correlations,  while  Von  Winkle  calls  them 
lateral  correlations.  If  it  is  assumed  that  the  wall  pressures 
are  normal  to  the  wall,  both  the  longitudinal  and  the  lateral 
(transverse)  pressure  correlations  are  of  the  lateral  type  as 
defined  by  de  Karman  and  Howarth,  since  the  components 
are  normal  to  the  direction  of  the  displacements.  These 
pressure  correlations  show  principally  how  the  fluctuating 
wall  pressure  is  related  along  the  wall  surface.  It  can  be 
expected  that  longitudinal  pressure  correlations  are  over 
larger  distances  than  the  lateral  pressure  correlations, 
since  the  main  flow  is  convected  in  the  longitudinal  (tunnel) 
direction. 

For  zero  displacement  distance,  the  nondimens ional 
correlation  of  shear  stress,  obtained  by  dividing  the  products 
of  the  rms  magnitudes,  has  been  mentioned  as  being  equal 
to  approximately  0. 45  over  most  of  the  pipe  radii: 


u  u 
r  X 


0.45 


The  correlations  of  this  report  are  autocorrelations. 
They  correlate  the  pressure  components  at  the  same  point 
and  in  the  same  radial  direction,  but  differing  in  the  time 
separation  or  lag  between  them: 


Mt)  = 
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liinT-*«  2i 


P{t)p{t^r)dt 


(28) 


(29) 


where  JR{r)  is  the  autocorrelation  between  the  fluctuating 
radial  pressures  displaced  by  the  time  lag  r 


T  is  the  time  interval  of  integration. 
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This  autocorrelation  of  pressures  is  similar  to  the  lateral 
correlation  between  radial  pressure  components  displaced 
in  the  direction  of  mean  flow:  the  displacement  distance  is 
the  product  of  the  time  lag  and  the  mean  flow  velocity. 


PERIODIC  FUNCTIONS  ANALYZED 

Since  most  water  pipe  flow  experimental  m.easure- 
ments  were  of  spectra,  a  brief  theoretical  discussion  of 
harmonic  analysis  follows;  it  is  a  generalized  Fourier 
series  approach  lo  the  study  of  many  functions.  A  Fourier 
series  is  an  infinite  series  expression  of  a  function  that 
repeats  itself  over  intervals  of  an  infinite  range  of  the  inde¬ 
pendent  variable.  The  function  is  expressed  as  an  infinite 
series  of  constants  multiplied  by  circular  functions  with 
periods  of  integer  multiples  of  the  period  of  the  original 
function.  The  constant  multipliers  are  obtained  from  the 
average  of  the  function  weighted  by  each  harmonic  circular 
function  over  the  fundamental  period,  because  of  the  im¬ 
portant  orthogonality  conditions  for  circular  functions.  The 
function  can  be  written  in  terms  of  the  exponential  form  of 
circular  or  sine  functions  with  both  positive  and  negative 
angular  velocities: 


At)  =  ^ 


n  =  0,  ±1,  ±2. . . . 
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where  T 
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fit) 
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is  the  fundamental  period 

is  the  periodic  function  expressed  as  an  infinite 
series  in  e 

-  in  ui  t 

is  the  complex  constant  of  the  series  in  e  '  i 
fit)  andFCn)  are  Fourier  transforms  of  each  other 

is  the  dc  coinponent  which  must  be  finite. 
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The  aistocorrelation  of  the  periodic  function  f[t)  is 
obtained  by  the  quadratic  average: 


1 

^(T)  *  f  f(mt^T)dt  - 
1 


n^-o) 


(31) 


where  |i^(n)j  is  the  absolute  magnitude  of  Fin) 
ji^(n)|®  is  the  power  spectrum  of /(t) 

Fir)  is  the  autocorrelation  for  lag  r. 

The  autocorrelation  and  the  po^^^er  spectrum  are 
Fourier  transforms  of  each  other: 
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(32) 


Both  the  Fourier  spectrum  and  the  power  spectrum  are 
discrete  functions  of  angular  velocity.  The  function  fit)  is 
a  periodic  function  of  time,  and  the  autocorrelation  is  a 
periodic  function  of  the  lag. 


APERIODIC  FUNCTIONS 

Aperiodic  functions  such  as  transients  are  analyzed 
by  the  Fourier  integral.  The  fundamental  period  of  the 
Fourier  series  is  allowed  to  approach  infinity,  and  an 
aperiodic  function  is  represented  by  the  Fourier  integral: 
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(33) 


where  /(t)  is  an  aperiodic  function  and  the  Fourier  trans¬ 
form  of  F((ji^ 

Fiui)  is  the  Fourier  transform  of  f(t) 

F(0)  is  required  to  be  finite,  che  condition  of 
convergence. 

For  aperiodic  functions,  the  autocorrelation  does  not  exist 
if  it  is  defined  as  above,  since  it  would  be  zero  when 
averaging  the  quadratic  over  an  infinite  interval.  The  power 
spectrum  does  not  exist  but  an  energy  spectrum  does  exist 
when  the  autocorrelation  is  defined  as  the  quadratic  sum 
without  dividing  by  the  period  to  form  an  average.  The 
function /(i)  is  aperiodic  by  definition,  and  the  Fourier 
transform  of  /(t)  is  a  continuous  function  of  frequency  or 
angular  velocity. 


RANDOM  FUNCTIONS 

The  water  flov/  noise  pressure  fluctuations  measured 
with  the  test  hydrophones  appear  to  be  random  functions  of 
time.  The  autocorrelation  function  can  be  assumed  to  exist 
and  be  related  to  the  power  density  spectrum.  The  auto¬ 
correlation  is  an  even  function  with  the  maximum  value  for 
zero  lag.  For  the  autocorrelation  function  to  be  fully 
»'ealizable,  the  power  spectrum  must  be  positive.  Both  the 
autocorrelation  fxmction  and  the  power  spectrum  are  real 
and  even  for  random  and  periodic  functions.  The  Fourier 
transforms  for  random  functions  can  be  written  as  cosine 
transforms: 


r 


r '4 


85 


T  ® 
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o 

where  G(f)  is  the  power  spectrum  based  on  only  positive 
frequencies 

B(r)  is  the  autocorrelation 

p(t)  is  the  wall  pressure  fluctuation. 

Parscval’s  theorem  is  obtained  by  setting  the  lag  equal  to 
zero.  This  corresponds  to  the  mean  square  value  of  the 
random  function  expressed  in  integrals  both  in  frequency 
and  lag  as  the  independent  variables: 
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where  p^  is  the  mean  square  value  of  the  random  function. 

The  numerical  integration  of  the  measured  power 
density  spectra  for  water  flow  noise  shows  that  they  are 
equal  to  the  broadband  magnitude  squared,  as  expected 
above.  This  shows  that  there  are  no  hidden  periodicies  in 
the  data. 


FUNCTIONS  SUMMARIZED 


Some  of  the  characteristics  of  the  transforms  are 
indicated  by  their  dimensions.  If  the  dimensions  of  f(tf 
are  assumed  to  be  of  the  dimensions  of  power,  the  auto¬ 
correlations  are  of  the  dimensions  of  pov/er.  The  power 
spectrum  is  in  the  terms  of  power  times  lag  time,  or  work. 
The  Fourier  integral  transform  of  a  function is  of  the 


86 


dimensions  of  f(t)  multiplied  by  time.  The  function  f{t) 
and  autocorrelations  must  be  aperiodic  if  their  transforms 
(which  are  in  the  dimensions  of  a  product  with  an  infinite 
time  interval)  are  to  be  finite.  Random  functions  do  not 
have  finite  Fourier  transforms,  but  their  autocorrelations 
are  aperiodic  and  their  power  spectra  are  continuous 
functions.  In  the  case  of  periodic  functions,  the  autocor¬ 
relation  and  the  power  spectrum  are  in  the  same  dimen¬ 
sions,  as  are  the  function  itself  and  its  Fourier  transform. 
Here,  the  function  and  the  autocorrelation  are  both  periodic, 
and  the  Fourier  spectrum  and  the  power  spectrum  both 
occur  at  the  same  discrete  frequencies.  The  nature  of  the 
transforms  of  various  function  types  is  summarized  in 
table  B-1. 


TABLE  B-1.  NATURE  OF  TRANSFORMS 
OF  VARIOUS  FUNCTION  TYPES 


1 

FUNCTION  TYPES 

1  t  1 

TRANSFORMS 

Periodic 

Transient 

Random 

Step 

1 

Impulse 

fit) 

periodic 

aperiodic 

continuous 

i 

continuing 

1 

discrete 

Fourier  spectrum 

discrete 

dc  component 
continuous 

infinite 

every\''/here 

infinite  at 
zero  dc 
component 
after  step 

continuous 

Power  spectrum 

discrete 

none  (energy 
spectrum) 

continuous 

discrete  at 
zero  only 

none 

Autocorrelation 

periodic 

1 

1 

none 

aperiodic 

continuous 

none 
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NUMERICAL  ANALYSIS 


AUTOCORRELATIONS 

In  the  water  pipe  flow  experiments,  photographs  were 
taken  of  oscilloscopic  tracings  of  hydrophone  signals,  to 
provide  time  records  of  sound  pressure  at  one  point  in  the 
pipe  for  numerical  evaluation  of  the  autocorrelations  and 
power  spectra.  The  form  of  correlation  coefficient  for 
numerical  data  is  the  sum  of  the  product  of  the  differences 
from  the  mean  values  for  the  tv/o  variables  divided  by  their 
magnitudes: 


where  x  and  y  are  two  variables  correlated. 


(36  > 


For  the  time  data,  the  record  is  divided  into  N  readings, 
x^t  r=l,  2,  3,  ...  equally  spaced  intervals,  r^,  with  the 
time  record,  T=Nt^  .  The  autocorrelation  between  the  values 
of  x^  at  two  points  on  the  record  is  obtained  for  the  lag  time 
between  them.  The  lag,  T,  is  taken  as  equal  to  ct  ,  with 
TO=0,  1,  2,  3,  /T  K<N,  The  autocorrelation  as  a  function 
of  the  number  of  lag  intervals,  w,  is  given  by. 
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If  5^.  and  are  essentially  constant,  the  autocorrelation 
equation  simpliftes  to  the  form: 


R 


A  -X  3 
m 


m  x^-x^ 


where  x 


1  ^ 

= 


1  ^ 
r=l 


The  power  spectra  were  calculated  by  an  equation 
that  does  not  use  the  autocorrelations  as  computed  above. 
The  number  of  time  data  in  the  photos  is  small,  even  though 
several  were  combined  in  the  analysis.  The  normalized 
power  spectra  v/ere  calculated  by  the  equation  equivalent  to 
the  numerical  Fourier  transform; 
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where  w  =  k(jt^-x^) 

f  =  — 

The  power  spectra  at  the  frequencies  corresponding  to  each 

n  were  smoothed  by  a  Hamming  method; 

Y~  ~  0.  23p  +0.54p  +0.  23p  ,,, 

^n-1  n+1 


Figures  5  and  6  show  the  autocorrelations  obtained  by 
the  numerical  method  from  the  photographs.  The  autocor~ 
relation  curves  give,  for  the  steady  flow,  the  Taylor  micro¬ 
scale  defined  by: 


1 


2 


The  integral  scale  is  obtained  from  the  integral: 


(40) 


=  f  Mt  (41) 

0 

With  uniform  moan  flow  velocity  u,  the  time  scale  in  the 
autocorrelations  is  equivalent  to  distances  equal  to  the  mean 
velocity  multiplied  by  the  time  values.  Approximate  average 
values  of  the  microscale  and  the  integral  scale  are  0.075  and 
0.11  inch,  respectively. 
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MICROSCALES 


Taylor  measured  the  microscales  of  turbulence 
which  he  himself  originally  defined.  His  longitudinal 
velocity  mj.croscale  is  a  function  of  flow  velocity,  being 
about  0.  315  inch  at  7.5  fps  flow.  The  lateral  velocity 
microscale  is  0.  707  of  the  longitudinal  for  isotropic  turbu¬ 
lence.  The  longitudinal,  or  tunnel-direction,  microscale 
has  an  average  value  of  0.15  inch  for  the  broadband  mea¬ 
surements  at  the  U.  S.  Navy  Underwater  Sound  Laboratory. 
The  lateral,  or  cross -channel,  microscale  has  an  average 
value  of  0.49  of  the  longitudinal  average. 


POWER  SPECTRUM 


Figures  9  and  10  show  spectra  numerically  evaluated 
from  the  photographs  by  the  procedures  mentioned.  Anal¬ 
ysis  of  random  functions  is  based  on  the  assumption  that  all 
samples  of  the  data  have  the  same  autocorrelation  under 
the  same  flow  conditions.  The  Wiener  theorem  states  that 
the  Fourier  transform  exists  if  the  autocorrelation  exists. 
The  power  density  spectrum  is  the  same  for  all  waveforms. 
The  power  spectra  were  also  calculated  by  the  numerical 
form  of  the  Fouriei  transform  of  the  autocorrelations: 
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h 

H-2  f. 

rm=l 


ft  cos2wvnf  T 

m  ''ll 


(42) 


m=l,  2,  3,  ...  ^  n=l,  2,  3,  ...  A 

where  P  is  the  power  spectrum 
n 

is  the  chosen  time  interval 
is  the  chosen  frequency  interval, 

1 


The  mean  squared  value  of  the  pressure  has  been  shown 
to  be  the  autocorrelation  for  zero  time  lag,  and  the  integrsil 
of  the  powei  spectrum: 
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(43) 


lim  T 
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Tl2 

J  p^Wt  = 

-T/2 


a 

/  GifW 

0 


where  G{f)  is  the  power  spectrum 
/  ifi  the  frequency 

T  is  the  time  period. 


ENERGY  SPECTRUM 

Graphs  in  this  report  show  the  square  root  of  the 
power  spectrum  plotted  against  wave  number.  Much  of  the 
theory  of  turbulence  has  been  concerned  with  the  frequency 
distribution  of  kinetic  energy.  Pressure  can  be  considered 
to  be  kinetic  energy  per  unit  volume,  at  least  in  dimensions. 
Thus,  the  square  root  of  the  power  spectrum  of  pressure  is 
of  the  same  dimensions  as  the  kinetic  energy  spectrum. 
Taylor  first  showed  the  transform  relationship  between  the 
kinetic  energy  spectra  and  the  correlations  of  the  velocity 
components.  Ifis  kinetic  energy  spectra  were  in  terms  of 
the  one -dimensional  wave  number.  The  wave  number  is 
equal  to  the  angular  velocity  divided  by  the  mean  velocity, 
and  has  the  dimensions  of  an  inverse  length.  Mathematical 
expressions  for  eddies  are  written  in  terms  of  three- 
dimensional  wave  numbers,  and  these  are  directly  related 
to  the  eddy  sizes.  Large  eddies  cause  fluctuations  at  low 
frequencies,  and  small  eddies  cause  high  frequency  fluctu¬ 
ations.  For  the  one -dimensional  wave  number,  the  inverse 
of  the  wave  number  cannot  be  directly  connected  with  eddy 
sizes. 


figure  B-12  indicates  the  general  form  of  the  kinetic 
energy  spectrum  plotted  against  wave  number  for  isotrop.c 
turbulence.  The  lower  wave-number  range  contains  the 
largest  eddies;  these  show  permanence  and  depend  upon  the 
condition  of  formation.  The  energy  spectrum  is  maximum 
in  the  range  of  energy-containing  eddies,  Tiie  reciprocal 
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H^SRGY  SPECTRUM 


LARC^ST  EDDIES  OF 
reRMANENT  CHARACTER 

EK!ERGY-CONTA?NING  INERTIAL 
EDDIES  SUBRANK 

CONST  »>1 


UNIVERSAL  EQUILIBRIUM 
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Figure  B-12,  Figurative  sketch  of  energy  spectrum 
levels,  indicating  relationships  (adapted  from  ref,  2), 
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of  the  wave  number  of  the  peak  energy  is  interpreted  as  the 
average  size  of  the  energy-containing  eddies.  In  the  highest 
range,  the  turbulence  becomes  statistically  steady,  with  the 
change  in  the  mean  values  becoming  small  with  time. 

In  this  Kolmogoroff  range  at  sufficiently  high  Reynolds 
numbers,  the  turbulence  is  independent  of  external  conditions 
and  is  in  statistical  equilibrium  determined  by  the  dissipation 
rate  and  the  kinematic  viscosity.  There  is  an  inertial  sub¬ 
range  for  very  high  Reynolds  numbers  where  the  energy 
spectrum  is  independent  of  the  kinematic  viscosity,  and 
depends  solely  on  the  dissipation.  The  predicted  slope  of 
the  energy  spectrum/wave  number  plot  in  the  equilibrium 
range  is  proportional  to  the  negative  5/3  power  of  the  wave 
n^’.mber.  Viscous  effects  become  important  as  the  wave 
number  increases. 
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